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ABSTRACT 


Polarized fluorescence from the cytoplasm of human lymphocytes 
was studied using a modified fluorescence spectrophotometer. This 
technique was developed by Cercek for measuring antigen/mitogen-induced 
changes in the average microviscosity of the cytoplasm and has been 
proposed as a diagnostic test for cancer (5). 

Experiments were performed with a Perkin-Elmer MPF-4 fluorescence 
spectrophotometer equipped with a thermostatted sample holder and a 
polarization accessory. 

The spectrophotometer method was found to be unreliable for 
accurately determining the polarization of fluorescence from the 
cytoplasm. Specifically, filtrate fluorescence was found to vary with 
the method of filtration. Syringe, suction and gravitational filtration 
techniques were tried. 

A Sates ape technique was then developed for determining 
the polarization from the cytoplasm. This method was found to be more 
consistent than the former but had limited application for antigen/ 
mitogen-stimulated systems. The method assumes that the polarization 
from the cytoplasm remains constant and that the rate of diffusion of 
fluorescein across the cell membrane is directly proportional to the 
intracellular concentration. Using the analytical method, the variation 
of the polarization of fluorescence from the cytoplasm with the pH, 


osmolarity and temperature of the supporting medium were found. 
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CHAPTER ONE 


INTRODUCTION 


1.1 CANCER DIAGNOSIS AND THE SCM TEST 


A test which reliably distinguishes cancer patients from those 
suffering from nonmalignant disorders, is simple to perform, and 
requires only a few milliliters of peripheral blood, would revolutionize 
cancer diagnosis. Such a test has been reported by L. Cercek, B. Cercek 
and C. 1. V. Franklin (5), and L.: Cercek, 8. Cercek antd.J.. V. Garrett 
(8), and is known as the structuredness of cytoplasmic matrix or SCM 
test. 

The SCM test was developed by the Cerceks (9-14) and employs the 
technique of fluorescence polarization. Using this technique, cells are 
irradiated with plane polarized light, and the polarization of 
fluorescence is simultaneously measured from fluorescein molecules in the 
cytoplasm. According to a classical theory by Einstein and Perrin (35), 
the polarization of fluorescence is related to the viscosity (n) by 
1/P « 1/n (a complete description is given in Appendix C). However, 
since the cytoplasm is a heterogeneous, polyphasic medium, the polariza- 
tion of fluorescence from the cytoplasm more accurately represents the 
average microviscosity or so-called "structuredness."" The extent of the 
cytoplasm for a typical mammalian cell is shown in Figure l. 

The application of the SCM test is based on the principle that 
stimulatory molecules iieraeeihe with receptors on the cell surface 


induce changes in the average microviscosity of the cytoplasm. Cercek 
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FIGURE 1. The extent of the cytoplasm for a typical mammalian cell. 


(From Brachet: Scientific American, 55 :205,1961) 
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et al. (8) found that they could distinguish between normal subjects and 
those suffering from chronic lymphatic leukemia on the basis of 
differences in the SCM following exposure of their lymphocytes to the 
mitogen phytohaemagglutinin (PHA). Normal subjects showed reduced SCM 
whereas cancer patients did not. Extending their studies to other 
malignancies, they demonstrated that lymphocytes from 41 patients known 
to have neoplastic disease of breast, bladder, bone, brain, ovary, 
uterus, larynx, oral cavity, lung or skin showed no decrease in SCM in 
response to PHA while lymphocytes from 70 of 71 healthy donors, and all 
of 17 patients with non-neoplastic diseases responded with a character- 
istic decrease (5). Field and Caspary (20) had previously shown that 
encephalitogenic factor (EF) and cancer basic protein (CaBP) could 
distinguish cancer patients from healthy volunteers using the macrophage 
electrophoretic mobility (MEM) test. Their assay also utilised small 
quantities of peripheral blood but was extremely difficult to perform 
and had the added disadvantage of failing to discriminate between 
malignancy, and neurodegenerative disorders. The Cerceks employed EF 
and CaBP in the SCM test and found the reciprocal effect to PHA: cancer 
patients lymphocytes showed a reduction in SCM but cells from normal | 
volunteers did not. This prompted the Cerceks to combine the responses 
to PHA and CaBP as a single parameter, the SCM response ratio (RRscp) - 
The RRogcy is the ratio of the polarization of fluorescence from the 
cytoplasm following exposure to CaBP to that following exposure to PHA 
(5). The ratio for donors without neoplastic disease ranges from 1.2 to 
1.7, whereas cancer patients have values between 0.6 and 0.9. Although 


either response appears to reliably detect malignancy, the Cerceks felt 
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the response ratio sharpened resolution. RRgcy values were used to 
monitor patients' lymphocytes after surgical removal of malignant 
tissues (6). Values were determined before and after operation. 
Surprisingly, 5 out of 6 patients showed a progressive increase in RRocwy 
which reached the healthy range within two weeks of operation. This 
suggests that the test could be a useful monitor of the adequacy of 
surgery. A further claim is the possible application of the test to 
screen for premalignant conditions (6). In 8 out of 12 patients with 
histologically declared benign growths of the breast, the RRgcy values 
were typical of malignant conditions. It was suggested that these 

might represent premalignant lesions. That the aberrant values from 
benign lesions questioned the interpretation placed on the results from 
malignant cases was not entertained. More recently, the Cerceks have 
demonstrated that decreased SCM can be elicited by stimulation with 
known tissue of the appropriate histological type (7). The specificity 
of response is illustrated by the demonstration that lymphocytes from 
patients with basal cell carcinoma responded to basal cell carcinoma 
tissue but not to squamous cell tissue. This suggests that the cells 
recognise tumour-specific antigens and broadens the diagnostic potential 
of the test enabling the type of malignancy to be determined. 

Although the Cerceks have not attempted a mechanistic explanation 
of the SCM phenomenon they clearly imply it has an immunological basis. 
By studying SCM values in single cells (15) with a microscope-fluorometer 
equipped with plane polarising filters, they have been able to show that 
15-23% of the lymphocytes from cancer patients respond to CaBP compared 


to 3-5% of the lymphocytes from healthy donors. The magnitude of the 
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response would appear to exclude participation of T (thymus-derived) 
lymphocytes through direct stimulation of membrane receptors for 
antigen, but could be compatible with the stimulation of cells having 
"passively" adsorbed antibody with antigen affinity. The difference in 
iO ccadiese response to PHA is less dramatic, normals having 45-54% of 
their lymphocytes showing decreased SCM compared to 36-45% of the 
lymphocytes from cancer patients. Since the overall response to PHA in 
cancer patients is negligible the degree of response per cell must be 
very small. Cancer patients frequently show normal PHA responses in 
the more conventional thymidine incorporation assay, suggesting that the 
SCM test is measuring a parameter unrelated to the blastogenic response. 

Despite the uncertainty surrounding the underlying principles, 
the diagnostic potential of the test is clear and warrants a stringent 
investigation. As of yet, no critical evaluation of the SCM test has 
been reported. It was therefore the purpose of this study to 
investigate the methodology of the SCM test and to standardize the 
procedure for possible clinical implementation. Specific objectives 
were: 

(i) To provide a standard procedure for measuring the 
polarization of fluorescence from the cytoplasm. 
(ii) To measure the average polarization of fluorescence from 

the cytoplasm of human lymphocytes and to obtain the normal variation. 

(iii) To establish the sensitivity of the average polarization 
of fluorescence with respect to variation in the temperature, pH, and 


osmolarity of the supporting medium. 
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CHAPTER TWO 


INSTRUMENTATION AND EXPERIMENTAL TECHNIQUE 


2.1 INSTRUMENTATION 


2.1.31. General Requirements 


The main difficulty in measuring polarized fluorescence from the 
cytoplasm of cells is due to Foérster transfer (Appendix C.2). Férster 
transfer is an energy exchange between molecules that depolarizes 
fluorescence and occurs when the concentration of fluorescent solute 
exceeds a certain critical value. In order to attain a measurable 
polarization from a system of cells, it is therefore important that the 
intracellular dye concentration does not exceed the critical amount. In 
general, this concentration is low (107+ M) and polarization measurement 
requires high optical sensitivity. 

The study of polarized fluorescence from biological systems is 
relatively recent and, as such, there are no instruments available 
specifically for this purpose. There are several existing devices that 
can be modified, however, and these include the fluorescence microscope, 
the cytofluorograf, and the fluorescence spectrophotometer. 

The fluorescence microscope enables the user to observe and 
quantitatively measure fluorescence from the cytoplasm of cells. This 
instrument can be modified for polarization study by adding two plane 
polarizing filters to the original optics. A schematic is shown in 


Figure 2. With the modified microscope, it is possible to measure the 
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polarization of fluorescence from single cells. Cercek (15) has 
recently done this with human lymphocytes using a modified Zeiss 
microscope-flourometer/01l. 

The disadvantage with this technique is that the microscope must 
be realigned for every cell that is examined. Consequently, the 
analysis of a large number of cells becomes tedious and impractical. 

The cytofluorograf, as an alternative, is designed to analyse 
“eels rapidly and automatically. A built-in laminar flow system provides 
typical sampling rates of 10° cells/minute (Bio/Physics Systems Inc.). 
The cytofluorograf is optically bi-functional and is primarily used for 
cell differentiation and counting. A diagram of the instrument is shown 
in Figure 2. To provide maximum sensitivity, a laser is used for 
excitation and, as each cell traverses the laser beam, it absorbs and 
scatters light relative to its morphology and stain uptake. At the 
choice of the user, the cytofluorograf simultaneously records either the 
fluorescent intensity at two different wavelengths (red and green), or 
the fluorescent intensity at one wavelength and forward angle scatter. 
In this way each cell is stored with two separate co-ordinates and can 
be later displayed on an oscilloscope. 

In spite of its potential usefulness for polarization studies, 
the cytofluorograf has limited optical sensitivity. This can be directly 
attributed to the rate of sampling and to the amount of signal obtainable 
from individual cells. For this reason, the cytofluorograf was not 


(1) 


chosen for instrumentation. 


lA group headed by Dr. A. Norman at the University of California 
is currently using a modified cytofluorograf to measure the polarization 
of fluorescence from cells. Direct correspondence (32) with this group 
confirms the sensitivity limitation of the instrument. 
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2.1.2 Fluorescence Spectrophotometer 

The most feasible method of measuring the polarization of 
fluorescence from a large number of cells is using a fluorescence 
spectrophotometer. Although this instrument is mainly used for 
luminescence spectroscopy, it is readily extensible to fluorescing 
biological systems. The spectrophotometer is designed so that for 
sufficiently dilute solutions (where secondary absorption is negligible) 
the recorded fluorescence intensity is directly proportional to sample 
concentration, 


de. I = Igp-K+Leeropec (2.1245 


where K is a constant of the instrument, Ip is the incident intensity, 

2 is the optical path length of the cuvette, and where e, op and C are 
the extinction coefficient, quantum efficiency, and concentration, 
respectively, of the sample. The quantities e« and $p depend on molecular 
structure and this topic is discussed in Appendix B.2. Most spectro- 
photometers employ right angle optics (emission recorded at right angles 
to the direction of excitation) and these can easily be modified for 
polarization study by inserting a plane polarizing filter before and 
after the sample holder. 

The main advantage in using the spectrophotometer is that the 
fluorescent intensity from a large number of cells (10°) can be measured 
simultaneously. This means that the average intracellular dye concentra- 
tion can be far lower than for single cell measurement which alleviates 
the problem of depolarization by Férster transfer. Additional advantages 
of the instrument include the optical versatility of separate excitation 


and emission monochromators and the facility of temperature control. 
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The difficulty with this method is that it introduces a background 
fluorescence comparable to the signal emitted by the cells. This is due 
to the amount of fluorescent dye that diffuses from the cells into the 
supporting medium. Consequently, considerable analysis must be performed 
to determine the polarization of fluorescence from the cytoplasm 
(Sections 2.4 and 2.5). 

Of the spectrophotometers commercially available, a Perkin-Elmer 
MPF-4 was chosen because of its sensitivity and overall versatility. 

This instrument was generously provided by Dr. John Charnock of the 
Department of Pharmacology. A diagram of the complete apparatus is shown 
in Figure 3. 

The source of excitation in the MPF-4 is a 150 watt continuous Xe 
arc. This light source is modulated at 56 HZ by an optical chopper to 
prevent 60 cycle interference during signal amplification. To enable 
monochromatic excitation and emission analysis, the instrument utilizes 
two precision Czerny-Turner monochromators. The reflection gratings 
comprising these are ruled vertically at 1200 lines/mm and are blazed 
at 300 nm. Each monochromator was calibrated with the 450.1 nm peak of 
the Xe source and both were accurate to within .5 nm. The entrance and 
exit slits for both excitation and Seaton monochromators are 
continuously adjustable and provide a spectral bandwidth from 0 to 20 nm. 
To eliminate scattered second order radiation, a selection of filters 
are provided preceding the emission monochromator. 

The photometric part of the instrument is designed so that a 
small portion of the excitation light is reflected by the beam splitter 


to a reference photomultiplier. The resulting electrical signal is then 
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FIGURE 4 . Sample compartment and polarization accessory for the 


fluorescence spectrophotometer. 
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13 
amplified and becomes part of a negative feedback circuit that controls 
the dynode voltage of the main photomultiplier. This feature enables the 
instrument to compensate for source intensity fluctuation. The main 
photomultiplier in the MPF-4 is an R446F wide range (200-900 nm) with a 
high red response. A separate amplifier controls the sensitivity of the 
instrument and amplifies the signal from the main photomultiplier. This 
unit also supplies the high voltage to the photomultipliers and contains 
the dynode voltage feedback circuitry. The output of the amplifier was 
displayed on a Perkin-Elmer 056 chart recorder. 

The temperature of the sample compartment was controlled with an 
external ethylene-glycol reservoir which was continuously pumped through 
the base of the sample holder. A Haake regulator was used to maintain 
the temperature of the external bath and cuvette temperature was 
monitored by a Fluke digital thermometer. 

The polarizing filters that were fitted in the sample compartment 
are manufactured by Polacoat and consist of a transparent quartz 
substrate coated with a 16 micron layer of a dichroic film. These 
filters provide a large working aperture (1'' diameter) and are designed 
to polarize efficiently through wide angles of incidence. An "'L" shaped 
holder positions the filters in the sample compartment and allows them 
to rotate through 90° (Figure 4). Each polarizer is calibrated in 30° 
increments from the pass direction which is indicated by a "0" on the 
outer edge. A restraining lever, which is fastened at 45° to the pass 
direction, allows the filters to transmit either vertically or 
horizontally polarized light in their extreme positions. The actual 


filter transmission is about 30% in the visible range. 
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14 
For all measurements done with the spectrophorometer, fused 
silica (Si0»), fluorescence free cuvettes were used. The cuvettes have 
a sample path length of 10 mm and are designed to hold a volume of 4 ml 


(12.5 mm x 12.5 mm x 45 mm outside dimensions). 


2.2 PREPARATION OF LYMPHOCYTES 


A normal adult has approximately 7 x 10° leukocytes per milliliter 
of blood. Since lymphocytes make up about 30% of all the white blood 
cells (Table I), this represents a nar possible yield of 2.1 x 10° 
lymphocytes per milliliter. The most common method of preparing 
lymphocytes is using a density-gradient separation (23). With this 
technique, the cells are separated by centrifugation according to their 
sedimentation rates. 

In general, the sedimentation rates of cells are determined by 
size rather than density and sedimentation is greatly enhanced if a cell 
type can be induced to aggregate. Solutions of Ficoll 400 are used for 
lymphocyte preparation because of its ability to promote erythrocyte 
aggregation. There are several Ficoll solutions commercially available 
and of these Ficoll-Paque (Pharmacia) was found to be the most effective. 
Ficoll-Paque is an aeuqous solution of density 1.077 + .001 g/cm? 
containing 5.7 g Ficoll and 9 g diatrizoate sodium in every 100 ml. The 
function of diatrizoate sodium is to provide the optimal density and 
osmolarity necessary for the efficient removal of other cells from the 
lymphocytes. 

To obtain the most homogeneous lymphocyte preparations, the 


following procedure was used: 
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AVERAGE WHITE BLOOD CELL PERCENTAGES 
FOR NORMAL ADULT (FROM GUYTON (22) ) 
% FROM TOTAL OF 
CELL TYPE 3 3 
7 x 10° WBC/mm 
POLYMORPHONUCLEAR NEUTROPHILS 
EOSINOPHILS 
BASOPHILS 
MONOCYTES 


LYMPHOCYTES 


IN ADDITION THERE ARE: 


Sie. 10° PLATELETS /mm> BLOOD 


5.4 x 10° RBC/mm> BLOOD 


TABLE I 
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Whole blood (10-20 ml) was collected in sterile vials and swirled 
over glass beads to catalyze clotting. Upon formation, the clots (and 
thus platelets) were removed and the remaining blood.was added to test 
tubes containing carbonyl iron (10 mg iron/ml blood). These tubes were 
agitated at 37° C for 10 minutes allowing granulocytes, monocytes, and 
macrophages to phagocytize the iron particles. The blood was then mixed 
with 3% plasma gelatin (75% blood, 25% gelatin) and left to stand at 37°C 
for approximately 30 minutes. During this time the plasma gel causes the 
erythrocytes to aggregate and these, along with the iron-saturated 
leukocytes, settle out under the force of gravity. 

The supernatant of the resulting suspension was gently layered 
onto Ficoll-Paque (3 ml Ficoll-Paque/3.5 ml blood in 1.5 cm diameter 
siliconized glass test tubes) and centrifuged at 400 g for 35 minutes at 
18° C. This process forces out the remaining erythrocytes and granulo- 
cytes and causes the lymphocytes to collect at the interface of the 
suspension system. 

The lymphocytes were then carefully extracted with a pasteur 
pipette and washed with Hank's balanced salt solution (3 ml lymphocytes/ 

7 ml Hank's). To complete the wash, the cells were recentrifuged at 300 g 
for 10 minutes and finally suspended at room temperature in Medium 199 
(Gibco). 

Cell counts were done with an autocytometer (Fischer) and were 
periodically checked with a haemocytometer. Cell viability was determined 
using the technique of fluorochromasia (37) and cell purity was determined 
by direct examination with a microscope. Typical yields of ./7 x 106 
lymphocytes/ml whole blood were obtained with this method and cell 


viability and purity were usually greater than 95%. 
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2.3 HYDROLYSIS KINETICS OF FDA BY LYMPHOCYTES 


2.3.1 Staining Technique 


In order to measure the polarization of fluorescence from a 
suspension of cells, it is necessary to stain the cytoplasm with a 
fluorescent dye. This can be done by immersing the cells in a solution 
containing a non-fluorescent flurogenic substrate. As the substrate 
diffuses into the cells, it is degraded by enzymes to yield a fluorescent 
end product. Rotman and Papermaster (37) have studied this process 
using various fluorescein esters of fatty acids. They found that the 
best substrate for staining was fluorescein diacetate (FDA). This ester 
is non-polar and rapidly diffuses across the cell membrane where it is 
enzymatically hydrolised to produce the fluorescent anion fluorescein 
(21). Because of its polarity, the outward diffusion of fluorescein is 
a much slower process than the inward flux of FDA. Consequently, the 
rate of production of fluorescein exceeds the rate of outward diffusion 
and the difference between these rates results in a net intracellular 
accumulation. 

Fluorescein is extensively used for staining because it has a 
high quantum efficiency (¢f = .80) and because it is not harmful to cells. 
This dye is particularly good for polarization studies because it does 
not complex to structures in the cell and is found uniformly distributed 
throughout the cytoplasm. Furthermore, fluorescein can be excited by 
visible light (470-475 nm) and irradiating cells in this range produces 


no damaging photolytic effects. 
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2.3.2 Michaelis-Menten Kinetics 


The enzymatic hydrolysis of FDA can be described by: 


FDA + EQS (FDA)E—F + E+ 2A (Zo 3.) 
2 


where (FDA)E is an intermediate enzyme-substrate complex and where F and 
A represent fluorescein and acetic acid respectively. For most enzyme 
reactions, the rate of formation of the enzyme-substrate complex is 
directly proportional to the free enzyme and substrate concentrations. 


Thus, the rate of change in the concentration of (FDA)E can be written: 


d{(FDA)E] _ 


ae K,-([E - (FDA)E]-[FDA] - (Kp + K3)+*[(FDA)E] (2.3.2) 


where [E - (FDA)E] is the free enzyme concentration. 
According to the Michaelis-Menten formulation (31), the reaction 


described by equation (2.3.1) quickly reaches a steady state. This 


d[ (FDA)E] 


at = 0, and the steady state concentration of (FDA)E 


occurs when 


is given by: 


> PFDA)? (ile Ug [EDA SLE] 
[(FDA)E] 5, = [FDA] + Ko + K3 = TFDA] + K, (253'53) 
Rj 


where Ky Rog is defined as the Michaelis-Menten constant. During 


the steady state, the rate of formation of fluorescein is described by: 


d[F] 
dt 


= K3+[(FDA)E] 5, (2.3.4) 
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and using equation (2.3.3): 


d{F] _ K3:[FDA]-[E] 


dt)  PeLrDAT +k ee 
Defining 
eC ein t -a[F] 2 
V Sales ee ° 
max " [ypA]>=™ dt. [E] eee, 


equation (2.3.5) can be modified to give: 


GE by 2 Vinax* [FDA] 


dt [FDA] + K, (2.3.7) 
i Kn 
ox: alF] = é det ihe, T C253.8) 
dt. Vmax [FDA] Vmax 


(Lineweaver-Burk modification (31)) 


Thus, for systems obeying Michaelis-Menten kinetics, the enzymatic 


hydrolysis of FDA can be completely described by the parameters K, and 


Vinx (equation 2.3.8), and these quantities can be obtained from a linear 
=] 
plot of ol against [FDA]~!. For cell systems, Ky and Vmax depend 


on the number of cells (enzyme) that are used for the hydrolysis reaction. 
However, having calculated K, and Vpax for a particular number of cells, 


the kinetics can be determined for any number by knowing the substrate/ 


cell ratio. 
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2.3.3 Experimental Determination of K, and Vay for Lymphocytes 

Because FDA is insoluble in water, stock solutions of 2.5 x 10°32 M 
FDA/acetone were prepared and stored at 4° C. Recrystallized FDA (Sigma) 
and 99% pure acetone ( Fisher ) were used for this purpose. Aqueous 
solutions (pH = 7.4) of FDA were made up by dissolving small volumes of 
the FDA/acetone solution in sterile phosphate buffered saline (PBS). 
Dilutions of approximately 1 in 10 produced minimal flocculation. 
Complete PBS (34, Table V) was used as a solvent for FDA/acetone 
because FDA is unstable in media with high amino acid content (Medium 
199) and immediately forms a highly fluorescent background. This effect 
makes it impossible to detect the fluorescein produced by the cells. A 
similar problem occurs when penicillin streptomycin is added to FDA/PBS 
and for this reason the antibiotic was deleted from the medium. 

Lymphocytes were prepared as in Section 2.2 and suspended in 
Medium 199 (Gibco) at concentrations of 5 x 10© cells/ml. Prior to 
incubation with FDA/PBS, the lymphocytes were kept in a water bath at 
Sie ue. 

For the kinetic study, the substrate concentration was varied 
from .25 uM to 10 uM FDA/PBS. Aliquots of .2 ml from the lymphocyte 
suspension were mixed with 2.8 ml of FDA/PBS and rapidly transferred in 
a cuvette to the sample holder of the spectrophotometer. The formation 
of fluorescein was recorded by irradiating the cells at 473 nm and 
monitoring the fluorescent intensity at 511 nm on the chart recorder. 
Measurements were taken for four minutes at 27° C. To determine the 
actual concentrations of fluorescein that were produced, a calibration 


curve was constructed using four fluorescein/PBS solutions of known 
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concentration. Since fluorescein is also insoluble in water, these 
solutions were prepared by first dissolving the fluorescein (Baker) in 
acetone and then dissolving small amounts of the acetone solution in PBS. 

The formation of fluorescein for five different substrate 
concentrations is shown in Figure 5. As can be seen from the diagram, 
the curves are not exactly linear and therefore Michaelis-Menten 
kinetics do not strictly apply. By fitting the experimental curves to 
various analytical functions, and using the "minimum chi-squared" 
criterion, it was found that the best fit was a quadratic polynomial. (2) 
The quadratic coefficient of these polynomials is most significant for 
the higher substrate concentrations (>10~© M FDA/PBS) and this indicates 
that the hydrolysis reaction is inhibited when the intracellular 
fluorescein concentration becomes too high. 

Rather than deleting the kinetics for substrate concentrations 


=] 
cats) against 


> 107° My, Km and Vmax were determined by plotting det ae 
[FDA]~!. A least squares fit to the data is shown in Figure 6 and from 


the X and Y intercepts K, = 5.65 x 107-© M and Vmax = 6.89 x 107!% moles/s. 
In Section 2.6.2, it is shown that after a four-minute incubation 

in FDA/PBS, and for a substrate/cell ratio of 7 x 107!° moles FDA/cell, 

about 50% of the total fluorescein Aodietd remains inside the cells. 

From Figures 5 and 6, this represents 2.02 x 1071! moles of fluorescein, 


or for a population of 10® lymphocytes, about 2.02 x 10°!” moles/cell. 


2For the analysis of the fluorescent intensity curves (Figure 5) 
a computer algorithm designed by Bevington (3) was used which provides a 
least squares fit to an arbitrary function. This program is listed in 
Fortran in Appendix D (Subroutine CURF). The experimental data were 
fitted to five different analytical functions including first, second, 
and third order polynomials, an exponential of the form I(t) = A(l-e7Bt), 
and the product of an exponential and a quadratic polynomial. 
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(moles x 10 ) 


FORMATION OF FLUORESCEIN 


0 2 4 
TIME (min) 
FIGURE 5 . Hydrolysis of FDA by lymphocytes. The reaction was 
15 


observed at 27° C for substrate/cell ratios of; (®) 1.4 x 10. ; 
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(sec /mole x 10°) 


d[F] 
dt 


| 


1 3 


ia (1 /mole x10) 


FIGURE 6 . Lineweaver-Burk plot for data appearing in Figure -5. 
Kn is the Michaelis-Menten constant and Vee is the maximum rate 


of hydrolysis of FDA . 
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Since the effective volume of a lymphocyte is 1.73 x 10> S51 this 
amount of fluorescein corresponds to an average intracellular concentra- 
Gidn ond, 177 xe 1077 M,C) 

As discussed in Appendix C.2, the polarization of fluorescence 
from fluorescein solutions begins to decrease for concentrations 
S 30%. Ot ME Thus, in order to attain an intracellular concentration 
below this value it is necessary to choose a substrate/cell ratio 


< 7 x 107!5 moles FDA/cell. 


2.3.4 Relative Fluorescence Yield 

Because of the relatively high cell concentrations that were used 
in these experiments, the fluorescence yield was studied to determine 
the effect of cell number. For this study, the cell concentration was 
varied from 1.67 x 10° to 1.33 x 10© lymphocytes/ml and the substrate/ 
cell ratio was maintained at 7 x 107!° moles FDA/cell. In each case the 
total amount of fluorescein that was produced from a four-minute 
incubation was appropriately scaled and compared to that for the 
3.33 x 10° cells/ml concentration (arbitrary standard). The results of 
the experiment are shown in Figure 7. From Figure 7, it appears that 
the fluorescent intensity decreases for concentrations > .67 x 10® cells/ 
ml. It should be pointed out, however, that the data include a 3-52 
scaling error and that the maximum loss in fluorescent intensity due to 


scattering was only 5%. 


3The average diameter of a lymphocyte is approximately 10 um and 
assuming the cells are spherical, this represents a volume of 5.24 x 
107-!3 1/cell. However, because the nucleus of a typical lymphocyte 
occupies 2/3 of the total volume, the actual volume of the cytoplasm is 
173° e120 8 eel. 
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RELATIVE FLUORESCENCE YIELD 
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LYMPHOCYTES / mi x 10° 


FIGURE 7 . Relative fluorescence yield. A cell concentration of 


SBE NS 2 10° lymphocytes/ml was chosen as an arbitrary standard. 
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The combined results of Figures 5, 6, and 7 indicate that the 
optimum substrate/cell ratio and cell concentration for polarization 


Study are 7 x 10-!5 ml/cell and 3.33 x 10° cell/ml respectively. 


2.4 MEASURING POLARIZED FLUORESCENCE 


2.4.1 Grating Polarization 


As described in Appendix C, the polarization of fluorescence is 
determined by exciting a See with plane polarized light and analyzing 
the fluorescent intensity bhawamlered /by a plane polarizing filter. 

More explicitly, the degree of polarization (equation (C.1.1)) is given 
by: 


ye ee ni 
hee eer: 


where I, and I, are the fluorescent intensities transmitted by a 
polarizer whose pass direction is aligned parallel and perpendicular, 
respectively, to the electric field vector (Eq) of the incident light. 

To obtain the maximum polarization using the spectrophotometer, 
excitation is done with vertically polarized light. Hence, the reference 
directions for measuring I, and I, are along the vertical and horizontal 
axes respectively. 

There is a minor difficulty using the spectrophotometer in that 
both I, and I, must pass through a grating monochromator before reaching 
the main photomultiplier (Figure 3). Consequently, Iy is significantly 
attenuated (relative to I,) by the vertical rulings of the grating. To 
compensate for the asymmetric attenuation Azumi and McGlynn (1) have 


devised a simple procedure for obtaining a grating correction factor. 
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Their method is to excite with horizontally polarized light and to 
record the ratio of the vertical and horizontal fluorescent intensities 
transmitted by an emission polarizer. Designating these components as 
Iyy and Igy respectively, the grating factor is gerined by G = Ipy/Ign- 
Because of the symmetry arising from horizontally polarized excitation 
(Appendix C) G should ideally be equal to 1 and any discrepancy is 
directly attributable to the grating. In practice G is always less than 
1 and depends on the wavelength setting of the emission monochromator and 
to a lesser extent on the emission bandwidth (slit width). For a 
particular value of G, the degree of polarization is correctly given by: 


Ty - Gly 
1g/PeeG=t, (22464) 


Pp = 

2.4.2 Experimental Procedure 

The general procedure for measuring polarized fluorescence from 
a suspension of lymphocytes was the same as for the kinetics study done 
in’ Section: 25313: 

Lymphocytes were prepared as in Section 2.2 and suspended in 
Medium 199 at a concentration of 5 x 10® cells/ml. Prior to incubation 
in FDA/PBS the cells were kept in a water bath at 37° C. 

To minimize depolarization by Férster transfer a substrate/cell 
ratio of 7 x 10-!° moles FDA/cell was used. Aliquots of .2 ml of the 
lymphocyte suspension were added to 2.8 ml of 2.5 x 10-© M FDA/PBS and 
transferred in a cuvette to the sample holder of the spectrophotometer. 
The cell suspension was irradiated at 473 nm with vertically polarized 
light (excitation polarizer in the "0" position as in Figure 4) and the 


fluorescent intensity components I, and I, were recorded at 511 nm on 
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28 
the chart recorder. Since I, and I, could not be recorded simultaneously 
the emission polarizer was alternated at regular intervals between the 
"oO" and "90" positions. Measurements were taken for approximately four 
minutes at 27° C. | 

A typical recorder printout is shown in Figure 8 where the data 
have been precorrected for grating attenuation. The solid segments of 
the I, and I, curves represent the actual experimental ie and the 
complete curves were obtained by fitting the data to quadratic polynomials 
(best least squares fit). The ce and If curves are hypothetical extra- 
cellular components to be described in Section 2.4.3. 

Figure 8 also shows that there is a small naturally occurring 
background fluorescence from the FDA/PBS solution. This component 
provided a convenient means of determining the grating correction factor 
and the background polarization. With the emission monochromator set at 
511 nm and using an emission bandwidth of 20 nm, G was calculated from 
the method of Azumi and McGlynn (1) to be .469 + .002. The value of G was 
checked periodically over the course of an experiment. 

From the I, and I, data the total polarization of fluorescence as 
a function of time can be calculated by: 
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2.4.3 Determination of the Polarization from the Cytoplasm 
Because of the diffusion of fluorescein into the supporting medium 
the total fluorescent intensity consists of an intracellular and extra- 


cellular component. By defining: 


2 Sa (2.4.2a) 
gud 40 feo ies als (2.4.2b) 


where e and i indicate extracellular and intracellular respectively, 


equation (C.1.1) can be expanded to give: 


Lge eevee ela pale vom tks 
=e i 4 ’ C2453) 
3 ie dy ye cl Pe i a ae 
This equation can be rewritten as: 
Diet lege Ph°Re (2.454) 


where P, is the polarization from the cytoplasm, Pp is the background 
Gh eee ; 
polarization, ud where Ry. = ssa such tthats2 Ro = Ke the is 
yn) aa ky n 
tempting to consider Ry, as the ratio of the intensity of the "n"th 
component to the total intensity. However, because It = Iy + 2I, 
(equation (C.1.12)) this is not exactly the case. To express the total 


polarization in terms of the actual component fractions, equation 


(C.1.10) can be expanded to give: 
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Equations (2.4.4) and (2.4.5) show that in order to calculate the 


polarization from the cytoplasm it is necessary to determine ie and ve 


2.4.4 Difficulties Obtaining P, by Filtration 

According to a method proposed by Cercek (12), P, can be 
experimentally determined by filtering the cell suspension at the end of 
a run and obtaining te and ie from the filtrate. By marking the time 
at which the cells are filtered on the chart paper and extrapolating Iy 


and I, to the time of filtration (te)s P. can be calculated from: 


Pict df ORR ENE SOE Ia 2.4.6 
ae (tpl oyit ty 2415) ce 
Although this technique appears straightforward, considerable 
dat ttoulty was encountered in applying it to the lymphocyte system. By 
trying various filtering techniques it was found that filtrate 
fluorescence Cie and La) is a function of the method of filtration. (4) 
Consequently, different values of P, are obtained corresponding to 
identical sets of Iy and I, curves. This inconsistency is the major 
difficulty in trying to determine P, by filtration. In addition to the 
filtering problem there is also significant error introduced by 
extrapolating I, and I,. Because of these difficulties, the experimental 


method for determining P, was abandoned. 


4Filtration was performed with 1" diameter Swinney filter holders 
(Canlab) containing 2 um polyvinylchloride filters (Millipore). Three 
separate techniques were tried including gravitational, syringe (3-5 ml), 
and suction filtration. For each method of filtration the filtrate 
fluorescence at t¢ = 4 min was distinctly different. The values of i 
and i were highest from the gravitational filtration followed by those 
from ee suction filtration and finally from those obtained by the 
syringe filtration. This indicated that the filtrate fluorescence 
depends partly on the force of filtration and also on chemical factors. 
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2.5 ANALYTICAL DETERMINATION OF Pe 


Since it was experimentally impossible to accurately determine 1¢ 
and If, an analytical method was developed. 

Extracellular fluorescein is the result of diffusion across the 
cell membrane and this process can be examined independently. Rotman 
and Papermaster (37) have studied fluorescein diffusion from various 
mammalian cells including human lymphocytes, and have found that the 
diffusion rate of intracellular fluorescein is proportional to the 


concentration remaining in the cell; 


= Kp: [F]4 C201) 


where [F] is the concentration of fluorescein, Kp is the diffusion 
constant, and where e and i refer to extracellular and intracellular 
respectively. This result was determined with a fluorescence microscope 
by observing the fluorescent intensity from single cells as a function 
of time. The results from a typical study are shown in Figure 9. 
Because of varying membrane properties, cells from the same population 
exhibit a great degree of heterogenéity with respect to diffusion 
constants. It is difficult to account for this in any analytical model 
and only an average diffusion constant can be assumed for a large number 
of cells. 

In Appendix C, it is shown that the total fluorescent intensity 
can be expressed by euqation (C.1.2) as It = Iy + 21,, and since I; is 


proportional to the total concentration of fluorescein 
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FIGURE 9. Release of intracellular fluorescein. A,B and C 
represent measurements taken from single mouse lymphoma cells 


(ML 388) in PBS at 23° C . (After Rotman and Papermaster (37)). 
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The total fluorescein concentration is the sum of the extra- 
cellular and intracellular concentrations and using [F]¢ = [Fle + [F]q 


equation (2.5.1) can be rewritten as: 


oree = Kp [Fla = Rp: (Fly — (Pe) 
oe a + Bp fFle = Kp [Flt (2.5.3) 
Following equation (2.5.2): 

Leah LS tao le Sali ie (2.5.4) 


and since equations (2.5.2) and (2.5.4) have the same constant of 
proportionality (assuming that the fluorescence efficiency of intra- 
cellular and extracellular fluorescein is the same), equation (2.5.3) 


can be expanded to give: 
d 
re Gly. t. 210) + Ry lps 219) = Kyi 4) 211) (235.55) 


This equation may be rearranged as: 


a4 a df f é 


which yields the two simultaneous equations: 


e .: 
diy 


az (Ut Kp*Iq = Kp'ly (2-506) 
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As stated in Section 2.3, the best least squares fit for both I, 


and I, was a quadratic polynomial and thus, for a lymphocyte systen, 


equations (2.5.6) and (2.5.7) can be solved for 15 and to. Since Iy 


and I, can be expressed generally as: 
Ty (t) = Ap + Ayt + Agt? C2. 


a particular solution for (2.5.6) is: 


(A, + 2Apt) 2Ad 


5.8) 


a5 


Pye = (Ay + Ayt + Aot2) - —————— + Fa 
TE6t) = Co 1 2t*) % ee (2. 2529) 
dif 
The complementary solution to the equation a: + Ky 1h = 0 is: 
cre (t) = Ae “pt (2.5.10) 
and the general solution to equation (2.5.6) is given by the sum of 
(2.5.9) and (2.5.10). At t = 0, I} (0) = Ao, and therefore 
A Ag Aq ihn fe Dhol Ry = Bo (2.5.11) 
Solving for A: 
A = Aj/Kp - 2A2/Kp (e519) 
and the general solution can then be written: 
A 2A = (A) + 2Aot) 2A 
e 1 _ 2Ag . 2, _ Ar + 2Agt) 2 
BeGt koe Kp ape? aGAgieh Anti ch Agt*) Kp Hees 
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A similar expression can be derived for Bett) - 
Equation (2.5.13) suggests an interesting method for calculating 


both P, and Kp simultaneously. By using equation (2.4.6), P,. is given by: 


tT, — 1p); -— 12) 
Prac Ey 4 seh ea Nuies dye 2s (2.5.14) 
Gato) (geet Gp Tp 


and an analytical expression for P.(t) can be obtained by simply providing 
a value for Kp: However, since the experimental conditions are well 


defined (pH = 7.4, isotonic media, constant temperature = 27° C), P,. 


dP 
should remain constant and therefore K should be chosen so that rah QO. 


dt 
This condition was impossible to achieve experimentally and Kp was chosen 
so that there was a minimum variation in P,. This was accomplished by 
calculating the average polarization <P,> for several values of Kp and 
finding the one with the minimum standard deviation. Because of the 
tedious nature of these calculations a computer program was drafted to 
determine the optimum value for <P,.>. The program was written in Fortran 
and is listed in Appendix D. For simplicity, a block diagram demonstrat- 
ing the main features is shown in Figure 10. The program is designed to 
calculate and list the average polarizations, <P>, and the standard 
deviations, op, for a range of diffusion constants supplied by the user. 
To eliminate random estimation of Kp, a subroutine CLCD1 was constructed 
which calculates Kp from a value of hig or a at a particular time. With 
this subprogram the user can initiate a search at a fictitious time of 
filtration and let the computer work its way through a range of possible 


filtrates try and i without. interruption. The optimum value of <P,> 


can then be chosen from the computer printout. 
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READS |, DATA 
FITS DATA TO QUADRATIC POLYNOMIAL 


READS |, DATA 
CORRECTS DATA FOR GRATING ATTENUATION 
FITS CORRECTED DATA TO QUADRATIC POLYNOMIAL 


READS I, AT t,> Al, , AND 
NUMBER OF DESIRED ITERATIONS 


FROM If 
FROM I¢ 


CALCULATES Kl 
CALCULATES Ki 


CALCULATES <P.> AND Op FOR 8 
SEPARATE DETERMINATIONS OF P. (t) 


PRINTS 
<P.> 


INCREMENTS 1, BY AI® 
PRINTS 8 SEPARATE 
VALUES OF P(t) FOR K, 


FIGURE 10 . Flow diagram of computer program for calculating Fe 5 
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2.6 EVALUATION OF THE ANALYTICAL MODEL 


2.6.1 Testing the Assumptions 


The analytical method for determining P, is based on two 


assumptions: 
d(F], 
(1) ie de Kp[F]l; (equation (2.5.1)) 
dP. 
C2) aes 0 (ee = constant) (2 sOreck:) 


To check the experimental validity of these assumptions two separate 
tests can be performed. 

Because i and Le are derived from separate quadratic poly- 
nomials, each function contains a distinct value for Kp (Ky and K)- 
However, according to equations (2.5.6) and (2.5.7), KR and ies should be 
identical, and therefore this criterion can be used to verify equation 
GP Ais Ti 

The second test that can be performed concerns the background 
polarization. Since the supporting medium is well defined (2.8 ml 2.5 uM 
FDA/PBS + .2 ml Medium 199 at 27° C), P}, should remain constant. Knowing 


Le and it this can be checked from the equation: 


i peMees ba 
i 
Pt) = (2.6.2) 
B re age ad 


The constancy of P, is a direct verification of equation (2.6.1). 
The "Ky" and a tests are a good means of evaluating a particular 


selection of <P.> and were therefore incorporated into the computer 
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program. The eh test proved to be very useful as large differences in 
the values of KS and ea were directly related to large magnitudes of op. 
In contrast, the "P," test was less significant because of the extremely 


low background polarization (Py = .015). 


2.6.2 A Typical Analysis 


To demonstrate the complete function of the computer program the 
numerical and graphical results from a typical run are shown in Figure 11 
and Table II. 

Figure 11 shows the variation a P, (equation (2.5.14)) for three 
values of Kp and demonstrates how the optimum value for <P.> is 
analytically determined. The total polarization (P) and the background 
polarization (P,) are also shown. 

From the numerical results in Table II, the optimum polarization 
was <P,> = .186 with a standard deviation of op = .002. On the basis 
that the per cent difference in Ky and ee was 9.9% and that Py = .014 + 
004 it can be concluded that .186 is a reasonable value for <P>. 

Also shown in Table II is the intracellular component ratio C4 
from equation (2.4.5). These data indicate that after a four-minute 
incubation about 50% of the total fluorescent intensity is from the 
cytoplasm of the cells. Since P, >> Py = 0, equations (2.4.4) and 


(2.4.5) can be approximated by: 
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. Graphical results from a typical run. 
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Equations (2.6.3) and (2.6.4) show that the time dependence of Cy and Ry 
is approximately the same as for the total polarization. The fact that 
Ci; gradually decreases can be explained by a steady-state intracellular 
concentration of fluorescein. When this occurs, i = constant, and 
since I; is continually increasing, i ve = C,; must decrease. From the 
data in Table II, the steady-state concentration is reached within the 


first minute of incubation. 
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CHAPTER THREE 


RESULTS AND DISCUSSION 


3.1 NORMAL RANGE FOR P. 


To establish the normal variation in Po for human lymphocytes, a 
total of 17 donors were tested. Data were compiled by calculating the 
average polarization ) and standard deviation Moe aerimene at? from 
approximately five runs for each donor. In all, 87 separate determina- 
tions of es were done and the results of the study are shown in Table III. 
The overall average polarization was calculated by weighting the 
individual averages according to their contribution to the total number 
of runs. This value is given at the bottom of the table as Pe EO ee 
.004. The overall standard deviation was calculated assuming a Normal 
distribution in the data and using an average standard deviation of .015. 
Havind determined the normal variation in P,, equation (C.3.27) 
was applied to determine the approximate range for the average micro- 
viscosity of the cytoplasm. Assuming that: 
(a) the principal polarization of fluorescein is 
Po = 1/2 (Appendix Gao) 

(b) the excited lifetime of fluorescein is 
t = 4.5 x 1079 seconds (48) 

(c) the molar volume of fluorescein is 


V5. = 000 mE (including solvation shell (35)) 


equation (C.3.27) reduces to (at 27° C): 
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NORMAL RANGE OF POLARIZED FLUORESCENCE 


FROM THE CYTOPLASM OF HUMAN LYMPHOCYTES 


experimental OF RUNS 


TABLE III 
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FIGURE 12 . Variation of the polarization of fluorescence from the 


cytoplasm of lymphocytes with the average microviscosity (27° C). 
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= 2.000 + .832 


3\|4 


al 
Py Ca24ed) 
(@ 


From a graph of this equation shown in Figure 12, ae value of Pa = .199 
corresponds to an average microviscosity of n= 30 centipoise (1 poise = 
1 dyne sec/cm*). This value is about 30 x more viscous than water at 
748 a Oe 

An average standard deviation of .015 represents about 8% of the 
value of hig From Figure 12, this represents a range in the average 
microviscosity between 25 and 35 centipoise. Although part of the 
standard deviation was due to fluctuations in tonicity and pH (Sections 
3.2.3 and 3.2.4), some of it was due to the analytical method for 
determining Pie SOE the data in Table III, the ee test showed that the 
majority of the % differences in KD and ss were in the 7-11% range. 
This corresponded to a variation in op between .002 and .010 and it was 
therefore estimated that the average contribution to Cevper imental from 


the analytical method was +.006. 


3.2 EXTERNAL FACTORS AFFECTING P. 


3.2.1 The pH and Osmolarity of the Supporting Medium 


In order to sustain mammalian cells "in vitro" such factors as 
the temperature, pH, tonicity and nutritional content of the supporting 
medium must be carefully controlled. Optimum cell growth can generally 
be attained by matching these parameters as closely as possible with 
those of the original cellular environment (extracellular fluid). The 


actual chemical compositions of the extracellular and intracellular fluids 
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47 
are quite different and are shown in Table IV. 
To maintain a pH in the vicinity of 7.4, the supporting medium is 
provided with a buffer system which consists of a weak acid in equilibrium 
with its salt. According to the Henderson-Hasselbach equation (31) the 


pH of a buffered medium can be described by: 


pil = pK, t log a Cre 2) 


where the acid dissociates as: 
16 Fema TS Wig (weak dissociation) 


Ande tnemsdilitadlssociaktes as: 


BAG oeeerecoeneee ca + A” (complete dissociation) 


As equation (3.2.2) shows, the pH of the medium will fluctuate around the 

value of pK, depending on the relative proportions of salt and acid. 
Medium 199 (Table VI) employs a bicarbonate/carbonic acid buffer 

which is the most widely used system for culture media. This system 


consists of a liquid-gas equilibrium described by the reactions: 


Ka 
+ = 
COn (ey + HO? 1002 a HCO, 


Giherel pk, = Gslvat 25° € (31)) and: 
NaHCO 3===2Na’ + HCO3". 


Applying equation (3.2.2) to these reactions: 


[NaHCO, ] 


pHe= Oke: TOPs TCO 
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CHEMICAL COMPOSITION OF EXTRACELLULAR 


AND INTRACELLULAR FLUIDS (FROM GUYTON (22)) 


EXTRACELLULAR | INTRACELLULAR 
COMPONENT FLUID FLUID 
pH=7.4 pH=7.0 
mEq/1 
mEq/1 
mEq/1 
mEq/1 


mEq/1 


HCO, mEq/1 


PHOSPHATES mEq/1 


SO, mEq/1 


GLUCOSE O - 20 mgm Z% 
AMINO ACIDS 200 mgm % 
CHOLES TEROL 


PHOSPHOLIPIDS . 2 oe ein, 
NEUTRAL FAT 


1 mEq = ‘hem Equivalent Weight 


1 mgm % = ioe gm/100 ml blood or plasma 


TABLE IV 
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49 
and can be controlled by adjusting the partial pressure of COo. A 
[NaHCO3]/[CO2] ratio of 20 to 1 provides a pH of 7.4. 
PBS (Table V) employs a phosphate/phosphoric acid buffer which is 
self-adjusting and is designed for short-term incubations. This system 


is described by the reactions: 


K 
a 
KH,PO, FST + Ht + HPO, 


(where PRI ©. 7 at 23" C Gl)) and: 
Na HPO, == 2Nat + HPO, ~~. 


From Table V, [Na jHP0O,]/[KH>P0,] = 4.59, and applying equation (3.2.2), 
pH = 7.4. The phosphate/phosphoric acid buffer has an advantage over the 
bicarbonate/carbonic acid system in that the value of pK, is closer to 
7.4. 

The simplest quantitative description of the tonicity of a solution 
is provided by the osmolarity (osmolality). The osmolarity is defined as 
the number of osmols per liter of solution where 1 osmol represents one 
mole of dissolved particles. (For example, 1 x 1073 moles of completely 
dissociated CaCl» represents an osmolarity of 3 milliosmols.) In 
calculating the osmilarity of concentrated solutions an osmotic 
coefficient must be introduced to allow for solute interaction (38). A 
154 mM solution of NaCl (isotonic saline) has an osmolarity of only .283 
osmols because the osmotic coefficient at this concentration is .92. 

Although all the components in the supporting medium contribute to 
the tonicity, the osmotic pressure is mainly controlled by NaCl. Fora 


solution of PBS the total osmolarity is .288 osmols while the osmolarity 


due to NaCl is .252 osmols. 
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CHEMICAL COMPOSITION OF 


PHOSPHATE BUFFERED SALINE 


~& (FROM PAUL@QG4 )) 


TABLE V 
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| | 1G MOLT TACT edi, 


* (CHOATE ii , 
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0.004 
G-00S 
Ovens, 
SOL 
(0008 


J 
‘ 


OsGeers 


Component 


AMINO ACIDS 

DL-Alpha-Alanine 

L-Arginine HC! 

DL-Aspartic acid 

L-Cysteine HC!I*H20 

L-Cystine 26.00 (2 HCi)(a) 
DL-Glutamic acid>H20 

L-Glutamine 

Glycine 

L-Histidine HCi*H20 
L-Hydroxyproline 

DL-Isoleucine 

DL-Leucine 

L-Lysine HCI 

OL-Methionine 

DOL-Phenylalanine 

L-Proline 

OL-Serine 

OL-Threonine 

DL-Tryptophan 

L-Tyrosine 57.88 (2 Na)(a) 
DL-Valine 

VITAMINS 

Ascorbic acid 

d-Biotin 

Calciferol 

Ca pantothenate 

Choline chloride 

Folic acid 

i-Inositol 

Menadione 

Niacin 

Niacinamide 

Para-aminobenzoic acid 

Pyridoxal HCL 

Pyridoxine HCI 

Riboflavin 
Thiamine HCl 
Vitamin A (acetate) (©) 


Remarks: 
(3)As supplied in powdered media 


(>)Omitted from 10X preparations and powdered media 


oe 


Component 


INORGANIC SALTS 


CaCl2 (anhyd.) 
Fe(NO3)3°9H20 
KCI 

KH2PO4 


MgSQq+7H20 
(anhyd.) (a) 


NaCl 
NaHCO3 (>) 
NaH2PO04*H20 


Na? HPO4+7H20 
(anhyd.) (a) 


OTHER COMPONENTS 
Adenine sulfate 


Adenosinetriphosphate (Disodium salt) 
Adenylic acid 

alpha tocopherol phosphate (sodium salt) 
Cholesterol 


Deoxyribose 
Glucose 
Glutathione 


Hypoxanthine (.354 Na salt) (a) 
Phenol red 

Ribose 

Sodium acetate 

Thymine 

Tween 80 At/as Powder Co. 

Uracil 

Xanthine (.344 Na salt) (a) ne. 


(c)Values established by the Tissue Culture Standards Committee. 


TABLE VI. The chemical composition of TC Medium 199 . 


( From GIBCO catalogue, 1976 ) 
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3.2.2 Physical State of H»0 in the Cytoplasm 

Two theories have been proposed to explain how the cell maintains 
the component distribution shown in Table IV. These are the membrane 
theory and the association-induction hypothesis (30). 

The membrane theory assumes that the membrane is the universal 
rate-limiting barrier to the traffic of water and all solutes between the 
cell and the external environment. Hypothetical pumps, located in the 
cell membrane are suggested to be responsible for excluding sugars and 
Nat or Cl~ ions from the cell. This theory also assumes that the bulk of 
cell water is normal and that the intracellular K* ions are in the free 
state. 

According to the association-induction hypothesis, a living cell 
represents a protein-aqueous fixed-charge system in which a high degree 
of molecular association exists among water, solutes and protein. This 
theory proposes that the component distribution in Table IV is maintained 
by metastable equilibrium which requires no direct energy expenditure. 

Kt ions are believed to be selectively adsorbed on protein anionic side 
chains and the bulk of cell water is presumed to exist as polarized 
multilayers. 

Recent NMR and diffusion studies (29) have provided strong 
evidence that the association-induction hypothesis is the most plausible 
of the two theories. The fact that the degree of polarization of 
fluorescence from the cytoplasm of lymphocytes is greater than 10X the 
polarization from the supporting medium is irrefutable evidence that the 
organization of the cytoplasm involves a high degree of molecular 


association. This result also suggests that the ratio of free to bound 
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water in the cytoplasm is an important factor in determining the average 


microviscosity and hence P,. 


53.2.3; The Effect of Osmolarity on P, 


It is an established fact that the osmotic pressure of the 
supporting medium controls the amount of intracellular water. To 
determine the effect of the osmolarity on P., the same procedure was 
followed as in ae 2.4 except that FDA/PBS solutions of varying 
osmolarity were used. Hypotonic and hypertonic solutions of PBS were 
prepared using NaCl concentrations of 4 g/1 (.159 osmols) and 16 g/1 
(.545 osmols) respectively. Solutions of isotonic PBS contained 8 g/l 
NaCl (.288 osmols). Approximately five determinations of P, were done 
for each osmolarity and the results are shown in Figure 13. 

In the hypertonic environment, where there is a net flow of water 
from the cells into the supporting medium, P, was found to increase to 
.223 + .010 at an osmolarity of .545 osmols. An opposite effect occurred 
for a hypotonic environment and P, decreased to .138 + .011. For isotonic 
solutions of PBS; Pe = <166 42 ).011. 

These results clearly indicate that the intracellular concentration 
of H20 is a critical factor in determining the polarization of fluor- 
escence from the cytoplasm. Similar data have been reported by Cercek (16). 

To determine the effect of small changes in the osmolarity a 
similar experiment was conducted by varying the concentration of CaClj in 
PBS. Hypotonic and hypertonic solutions were prepared using CaClo 
concentrations of 0 and 2 x 10-3 M respectively. Isotonic solutions of 
PBS contained ~ 1 x 107° M CaClp. Again, approximately five determinations 


of P, were done for each osmolarity and the results are shown in Figure 14. 
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FIGURE 13 . Variation of the polarization of fluorescence from 


the cytoplasm with the osmolarity of the supporting medium. The 


osmolarity of isotonic PBS is .288 osmols . 
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Because the variation in the concentration of CaCly did not significantly 
affect the total number of osmols in PBS (2 x 10-3 M CaClo corresponds to 
a change in the osmolarity of .0027), the units of the x axis of Figure 14 
are shown as concentration. 

For a CaCl2 concentration of 2 x 107° M, P. was found to be 
~-192 + .012 while for zero concentration P, = .170 + .011. The value of 
P. corresponding to isotonic PBS was P, = .184 + .011. Although these 
results indicate that P, is sensitive to extremely small changes in the 
osmolarity it should be pointed out that there is considerable overlap in 
P. for all three osmolarities. 

It was interesting to note in these experiments that the percentage 
difference in Kh and KH for the individual determinations of P, was in 
the same 7-11% range as for the normal study (Table III). This indicated 
that the exchange of H20 between the cells and the supporting medium 
occurred within the first minute of incubation in PBS. If the exchange 
occurred over a period of several minutes, it would be expected that P, 
would not remain constant (violation of equation (2.6.1)) and that there 
would be unusually large differences in Ky and Kp. It should be 
emphasized that the analytic method for determining P, depends on the 
validity of equation (2.6.1). 

During an experiment, it is unlikely that the osmolarity of the 
supporting medium changes appreciably. It can therefore be concluded 
that fluctuations in the osmolarity can only account for a small 


. . It is possible, however, that changes in 
percentage of Coxperimental Pp : ’ & 


the osmolarity on a day to day: basis could account for a significant 


variation in the magnitude of P.. 
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FIGURE 14 Variation of the polarization of fluorescence from 


the cytoplasm with the concentration of CaCl, in the supporting 


medium. The concentration of CaCl, in isotonic PBS is 
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3.2.4 The Effect of pH on P, 

Since changes in the tonicity of the supporting medium can 
significantly affect P., it was expected that variation in the pH would 
produce similar results. To study this effect, the same procedure was 
followed as in Chapter 2.4 except that the cells were immersed in 
isotonic PBS of varying pH. 

The effect of increased pH could not be determined because FDA is 
unstable under these conditions (pH > 7.4) and decomposes into fluorescein. 

The effect of decreased pH was studied by using three solutions of 
PBS with pH's of 7.4, 6.5, and 5.5. The acidic solutions were prepared 
by adding .1 N HCl dropwise to isotonic PBS until the pH stabilized at the 
appropriate values. These were read on a Radiometer pH meter/28 and were 
accurate to within +.05. Approximately seven runs were done at each pH 
and the results of the experiment are shown in Figure 15. 

As the pH of the medium decreases, P, decreases from .186 + .014 
BG pH) = 7.4 to .L75:* .0L0 at pH = 6.5 and? finally te: 141 = .01L at: pH 
= 5.5. Figure 15 also indicates that P, does not increase significantly 
for pH > 7.4. 

Although the pH effect on P, appears to be similar to the hypotonic 
effect described in Section 3.2.3, the exact mechanism is more complicated. 
Legge and Shortman (39) have previously reported that erythrocytes and 
thymic lymphocytes swell when immersed in acidic media. These authors 
proposed a mechanism for this process whereby HCO3 and Cl™ ions first 
enter the cell according to Gibbs-Donnan equilibria and are then followed 
by an influx of water in response to a hypertonic intracellular fluid. 


The data in Figure 15 are compatible with this mechanism and demonstrate 
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FIGURE 15 . Variation of the polarization of fluorescence from 


the cytoplasm with the pH of the supporting medium. 
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again the importance of intracellular H20 in determining the average 
microviscosity and hence P,.. 

The percentage difference in ms and KD for this experiment was 
also in the 7-11% range which indicated that the analytical method for 
determining P, is valid for acidic media. 

Because the pH is the most difficult parameter to control 
experimentally it is highly probable that most of Sexperimental is due to 
fluctuations in the pH of the supporting medium. The pH effect could also 
account for day to day variation in the magnitude of P.. 

The effect of pH and osmolarity on P, can be summarized by 


equation (C.3.27) as: 


bibs 


where n = f (pH, osmolarity) at constant temperature. 


3.2.5 Effect of Temperature on Po 


Prior to this section, no explanation has been given as to why P, 
was measured at 27° C. At temperatures >28° C, the kinetics for the 
hydrolysis of FDA are extremely unpredictable and sporadic increases in 
the fluorescent intensity make it impossible to accurately fit I, and 
I, data to analytical curves. At temperatures <28° C, the kinetics are 
uniform but the efficiency of enzymatic hydrolysis (and hence the 
magnitude of the I, and I, components) decreases with decreasing 
temperature. Therefore, an intermediate temperature of 27° C was chosen. 

To simulate "in vivo" metabolism, the lymphocytes were kept at 
37° C prior to their immersion in FDA/PBS. It should be mentioned that 


prolonged incubation (>is hr) at 37° C caused the cell médium to turn 
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alkaline and resulted in non-uniform kinetics. 
From the best data that could be obtained, it was found that the 


temperature dependence of P. follows that predicted by equation (C.3.27), 


Pee. 5 ue 


where n = f (T) decreases slowly with increasing temperature. For a 
temperature = 30° C, P, was determined to be 7% less than the average of 
.199 while at room temperature (23° C), P, was consistently 5% higher than 
oo. 

Because the temperature was controlled to within +.6° C (including 
a temperature gradient in the cuvette) it can be concluded that the 


experimental variation in temperature does not account for a significant 


part of Oexperimental- 
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CHAPTER FOUR 


SUMMARY AND CONCLUSION 


In 1974, Cercek et al. (5) reported that they could reliably 
distinguish cancer patients from those suffering from non-malignant 
disorders with the SCM test. The SCM test is performed on lymphocytes 
and employs the principle that stimulatory molecules, interacting with 
receptors on the cell membrane, induce changes in the average micro- 
viscosity of the cytoplasm. Such changes can be detected indirectly by 
measuring changes in the average polarization of fluorescence from 
fluorescein molecules in the cytoplasm (P,). 

The original objective in this study was to investigate 
applications of the SCM test but sufficient difficulty was encountered 
establishing a well defined control that the emphasis has been on the 
methodology. Although the SCM test was originally performed with a 
fluorescence spectrophotometer, this study suggests that such a method 
is unreliable for detecting changes in the average polarization of 
fluorescence from the cytoplasm. The principal evidence supporting this 
claim is that the filtrate fluorescence le and i) was found to vary 
with the method of filtration (Chapter 2.4.4). Since the Cerceks' 
procedure for determining P. depends on the extracellular or filtrate 
fluorescence, it is clear that the method is inconsistent. In addition 
to the filtration problem, fluctuating pH and data extrapolation are also 
significant sources of error in determining P,. 

In an attempt to circumvent the filtration difficulties, a semi 


analytical method was developed for determining P,. (Chapter 2.5). It 
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should be emphasized that the method is limited to systems where P, 
remains constant and where the rate of diffusion of fluorescein across 
the cell membrane is directly proportional to the intracellular 
concentration. For antigen/mitogen stimulated systems it is questionable 
as to whether these constraints are applicable. Applying the method to 
unstimulated lymphocytes, the average polarization of fluorescence from 
the cytoplasm (isotonic PBS at 27° C) was determined to be Pe = .199 + 
-004. This value is in accord with that determined by the Cerceks and 
corresponds (Perrin equation) to an average microviscosity of n = 30 cp. 
The analytical method can also be used to demonstrate that P, is extremely 
sensitive to variation in the pH or osmolarity of the supporting medium 
(Chapter 3). This indicates that the intracellular water content is a 
key factor in determining the average polarization of fluorescence from 
the cytoplasm. According to the Cerceks' data, the value of P, decreases 
by approximately 20% as the result of antigen or mitogen stimulation. It 
is interesting to note that the same decrease can also be achieved by 
suspending the cells in a medium with a pH = 5.5 or with an osmolarity of 
about .2 osmols. These facts suggest there are two separate mechanisms 
by which P, can decrease. First, by changes in the intracellular water 
content as the result of osmosis, pov second, by molecular rearrangement 
or structural transformation within the cytoplasm following antigen/ 
mitogen stimulation. 

It should be pointed out that this study does not directly refute 
the antigen/mitogen-induced effects reported by the Cerceks. To properly 
evaluate such phenomena, a different technique than the cell suspension 


method (spectrophorometer) must be used. Although the limitations of 
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alternative instrumentation were discussed in Chapter 2.1, it appears 
that a highly sensitive fluorimeter with a built-in flow system is the 
only feasible solution. Because of the obvious potential associated 


with the SCM test such systems must be investigated in the near future. 
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APPENDIX A 


POLARIZED LIGHT 


Light, or electromagnetic radiation, can be described as the 
superposition of plane waves.! These consist of sinusoidally varying 
electric and magnetic fields aligned perpendicular to each other and to 
the direction of propagation (Figure 16). 

To facilitate the discussion and pictorial representation of 
these waves it is customary to refer to only one of the two fields and 
usually the electric. Knowing either of the fields the other may be 


deduced from Maxwell's equations. 


A.1l Pure Polarized Light 


For a single monochromatic plane wave travelling in the positive 


z direction with respect to an xyz Cartesian coordinate system, the 


las derived from Maxwell's equations, the three-dimensional wave 
equations governing the behavior of the electric and magnetic fields in 
linear isotropic nonconducting media are: 


ae 2% 
v2 ot? v2 9t? 
bg > > 
With ve= c/s in = Ce) “3 Ds chsh = un 


where c is the wave speed in a vacuum, v is the phase velocity, and where 
n, € and uw are the index of refraction, dielectric constant, and magnetic 


permeability of the medium respectively. 


Analytic plane wave solutions to these equations (25) are: 


>> : >> 
i= E eifwt - eer = Bo ai(wt - k r) 
with k.= ra = w/v, where w is the frequency of oscillation of the fields. 
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PLANE WAVE SOLUTION V4 
TO THE WAVE EQUATION 


PLANE POLARIZED LIGHT 


FIGURE 16 . Spatial and temporal cross sections for plane 


polarized light. 
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physical electric field is given by: 


i(wt - kz)y 


E(x,y,z,t) = Re{E,e (A.1.1) 


where Re{ } and Im{ } represent the real and imaginary parts respectively 
of a complex number, and where Eo is a vector in the xy plane that may be 
complex; 


ee, E, = (E.,E.,0) (ar pe?) 


where E =A goo and E =A ely, 
te x. iv8y 


The individual components of the total electric field are therefore: 


E.. = A,cos (wt - kz+ b,.) (Ac Bsa) 
EY = ApeOs (OE - kz + *,) CAL 1636) 
E. = 0 CA. Le 3c) 


Different types of pure polarized light are distinguished by the 


relative sizes of AY and A, and on the phase difference > defined by: 


= 0 (A.1.4) 


In general the quantities BE ene and » assume arbitrary value and 
the light is elliptically polarized. From a fixed point on the z axis 
the electric field vector traces out an ellipse in the xy plane with an 
angular speed w. The ellipse is contained in a rectangle whose sides are 
parallel to the coordinate axes and whose dimensions are 2A, byg2Ate The 
eccentricity of the ellipse e« depends on the magnitude of » while the 
helicity, or sense of rotation about z, depends on the sign. LE to> 30 
then the sense of rotation is clockwise (negative helicity) and the light 


is defined to be right elliptically polarized (REP). Conversely, for 
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o < 0 the rotation is counterclockwise (positive helicity) and the light 
is left elliptically polarized (LEP). For fixed t the electric vector 
traces out a helix along an elliptical cylinder parallel to the z axis 
with a cross section identical to the ellipse referenced above for fixed 
z. The sense of rotation of this helix about the z axis is opposite to 
that of the electric vector for fixed z. These features are shown in 
Figure 17. 

Plane polarized light occurs for » = 0 or » = 7 and the electric 


field components from (A.1.3) are given by: 


€A.1°5) 


where EY is positive for > = 0 and negative for > = 7. 

For either of these cases the eccentricity of the ellipse describing the 
electric field is a maximum (ce = 1) and for fixed z the electric field 
undergoes simple harmonic motion along the diagonal of the rectangle 
(Figure 16). At a given time, the electric field undergoes simple 
harmonic motion in a plane containing the diagonal of the rectangle and 
running parallel to the z axis. 

When ¢ = + 5 and AL ¢ A the eccentricity of the ellipse is a 
minimum (0 < ¢€ < 1) and the major axis of the ellipse lies along the x 
axis. 

For ¢ = + 5 and A = A_ the eccentricity is zero and the light is 
right or left circularly polarized (RCP or LCP). In these cases the 


electric field components from (A.1.3) are: 
E = AS cos(wt - kz + $,,) (A.1.6a) 


Eo = AD sin(wt - kz + o,) (A.1.6b) 
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RIGHT ELLIPTICALLY 
POLARIZED 


LEFT CIRCULARLY 
POLARIZED 


FIGURE 17. Spatial and temporal cross sections for elliptically 


and circularly polarized light. 
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where He is positive for $¢ = - 5 and negative for 9 = + = This is 
shown in Figure 17. 
Although elliptically polarized light may be described as the 
superposition of two orthogonal plane polarized components an equally 


valid depiction is the superposition of left and right circularly 


polarized components (28). Defining 6: 
8 = wt - kz CAL Ey) 


the LCP components from (A.1.6) can be written: 


E, = A, cos(@ + $,) = Re{A, e 9 e+} (A.1.8a) 
B, = A, sin(@ + 6,) = Ima, otal Pass (A.1.8b) 
and similarly for the RCP components: 
-i60 id¢_ 
E.. = A_cos(@ - ¢_) = Re{fA_e e } (A.1.9a) 
B, = - A_sin(® - ¢_) = ImfA_ Gon tent) (A.1.9b) 


With equations (A.1.8) and (A.1.9) the total analytical electric field 


is written: 


B= 2, + iB, = quent es Ae ?= et (A.1.10) 
By defining a and 8 such that: 
a (A.1.11a) 
a aa (AL .11b) 


equation (A.1.10) can be written: 
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Bee oe: rey on ee al?) (A.1.12) 


Equation (A.1.12) completely defines the general case of elliptically 


polarized light. For 6 = a/2, 


E = aa +A) (A103) 


and the magnitude A, + A_ (which is the physical electric field) is the 
major axis of the ellipse making an angle 8/2 with the x axis. For 


6 = a/2.+ n/2, 
ee tit wore 5 (A.1.14) 


and the magnitude A, - A_ is the minor axis making an angle 6/2 + 1/2 


with the x axis. These features are shown in Figure 18. 
Finally, the essential parameters from both the plane polarized 


and circularly polarized superpositions are related by: 


sin 2n = sin 2y sin $6 (Ac 1 .a5a) 


tan 8 = tan 2x cos ¢ (A. 14. 15b) 


where tan 8 is positive for REP and negative for LEP, and where: 


1 A 
a tan. G) CAG 6) 
x 
n = tan “(|A, - A_|/(A, + A_)) (A.1.17) 


Ae: spare tally Polarized Light 


In Section A.1 pure polarized states were described for single 


monochromatic waves. Since light can be described as the superposition of 
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FIGURE 18 . Elliptically polarized light as the superposition of 


left and right circularly polarized light. (After Klein (28) ). 
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such waves, the polarization characteristics of a macroscopic beam of 
monochromatic light can be defined in terms of the time averages of two 
perpendicular field components and the mutual coherence of these 
components. 

However, in contrast with pure polarized states, where there is a 
well defined phase difference between two perpendicular field components, 
the superposition of many waves can produce a phase difference that 
changes rapidly over the period of observation and remains undefined. 
When this occurs, the time average of ae two electric field components 
(in a plane perpendicular to the direction of propagation) will be the 
same, and the light is said to be unpolarized. 

If unpolarized light combines with pure polarized light the 
resulting mixture is partially polarized. To characterize partially 
polarized light the degree of polarization is defined as the ratio of 
the intensity of the polarized part to the total intensity (26). 

To physically determine the state of pure polarized light or the 
degree of polarization in the case of partially polarized light, a plane 
polarizing filter and a phase shifter are ordinarily needed. This 
procedure is described in Section A.3. 

A plane polarizing filter converts unpolarized light to plane 
polarized by transmitting only the electric field components aligned in 
a particular direction. This is referred to as the pass direction. 
Polarizers are usually made with a dichroic material or less commonly 
with birefringent (doubly refracting) crystals (40). 

A phase shifter, or compensator, acts on a normally incident beam 
by retarding the phase of all electric field components aligned 


perpendicular to a fixed axis relative to the field components aligned 
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parallel to the axis. The direction along which the field components are 
retarded is called the slow axis while the fixed axis is called the fast 
axis. A compensator imposing a phase shift of 1/2 is called a quarter- 
wave plate. By positioning the fast axis of a quarter-wave plate along 
the direction of the major axis of elliptically polarized light the 

light is converted to plane polarized. This follows from the fact that 
the resultant phase difference between field components parallel and 
perpendicular to the major axis will be O or tm depending on the helicity 
of the incident light. Compensators are made from uniaxial birefringent 


crystals (40). 


J 8h" IN Coherency Matrix 


Recalling equation (A.1.1), the analytic electric field components 


at a particular point in space are given by: 


Egebial eifwt + 4) ChE Sia) 


ll 
I 


a ei(ut + $y) (A. 3.1b) 


By placing a compensator in the light with its fast axis along the 


x axis the field components become: 


E =A ei (ut a = (A.3.2a) 
x x 


=F ei (ut + Gr =. ) i ie et 7 (A. 3.2b) 
y », Mf 
where 6 is the amount that oy is retarded relative to d, 

If the resulting beam then passes through a polarizer with its 


pass direction inclined at an angle a to the x axis, the emerging signal 
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will be the sum of the projections of Ee and EY along a: 


E(a;6) = Ecos a + ants EY sin a (Re Sa3) 
The intensity? can then be expressed: 


I(a,6) = < E(a,6) E*(a,S) > CK23).8)) 


where the notation < > and * represent time average and complex conjugate 
respectively. 


Expanding (A.3.4) in terms of (A.3.3) the intensity becomes: 


I(a,6) = < E EX > cos¢a + < E EX > sin2a + 
p.ame« Vian; 


(</E EX > Behe + < EXE > ones sin a cos a (A.325) 
x y x y 


Following a similar derivation by Klein (28), the four time averages 


appearing in (A.3.5) may be described as the elements of a 2 x 2 matrix 


2Although not technically correct, the intensity is described: 


1 = |E,|? 


which is an abbreviated form of the mean irradiance. The irradiance, or 
electromagnetic flux density, associated with the plane waves described 
above is given by (25): 

> > > 

S=ExH 
where $ is Poynting's vector. The mean irradiance can be written: 


Senne Se =F ble gl4 


where the factor of % evolves from averaging a cos* function of time. 


Since the derivation of the degree of polarization involves the 
ratio of two irradiances, the expression for I as given above may be used. 
I is only used to simplify the derivation. More physically correct: 

> 
2 
I « |E,| 
(MKS units have been used in the discussion above.) 
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called the coherency matrix: 


= * = Tt as 
ee : re ‘ Cay i BE” 
C = 
Babe te (A.3.6) 
C = < EXE > C = <1 Be > 
yx x y yay yao, 
Since ce = ae these quantities may be expressed as: 
Cau WiC | ea (A.3.7a) 
xy xy e e 
Ce Slay.) 367? (A. 3.7b) 
yx xy 
Equation (A.3.5) can then be simplified: 
= 2 fine 
T(a,6) Cx cos‘a + Coy sin“a 
+ 2\c,,| sin a cos a cos(6é - $) (A.3.8) 
When (A.3.8) is averaged over a: 
< 1(a;6) cee Dye Cx at C7 = 1/2 I, (A.3.9) 


where I. is the total intensity of the beam. The total intensity can be 


written: 


I. = Cs te Cy =o Tr .¢ CASS 20) 


where Tr indicates the trace or sum of the diagonal elements of a matrix. 
Using equations (A.3.1), the coherency matrix for arbitrary elliptically 


polarized light is written: 
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< A@ > e i¢ <5 AGV AVS 
x x y 
ae CH. 381) 
el? SARIN, He Kes 
x y Vy 


and using equation (A.1.16) with I, = < AZ >+< AS >: 


cos*x esi? sin x cos x 
ai (A'.3.12) 
ig 


e * sin x cos x sin*y 


Examining C shows that det C = 0 for all pure polarized states where 


det C=C. ¢.. = Cur. 
Xe yy: vx xy 


For unpolarized light there is no coherence between the x and y 


C 


components and therefore hes = Cox = 0. Also, since the time averages 
of any two perpendicular components must be equal, < A, >=< ay >, and 


the coherency matrix for unpolarized light is written: 
i 


C= + (K53.13) 


These matrices are used to calculate the degree of polarization of 


partially polarized light. 


A.4 The Degree of Polarization 


An arbitrary coherency matrix may be written as the sum of a 


polarized and an unpolarized matrix: 


Gar + = c' + cP (A.4.1) 
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where u and p indicate unpolarized and polarized light respectively. 


C may also be expressed: 


C =atb =d 
xx xy 
C= CA a2) 
Cc = d* Ce td 
yx Wy) 
From section A.3, det cP = 0 so that: 
be - d¥d = 0 CAS 4593) 
Substituting b = Cs - a, and c= ae - a from (A.4.2) into equation 
(A.4.3) yields: 
aeve al (Cl) Ome. Cus = Oe SO (A.4.4) 
xx P yy AGG AV) Seay, YH 
or restating (A.4.4) in terms of the trace and determinant of C: 
a> - a Te. C + det C = 0 C&C 4.5) 
Following Klein (28), the only meaningful solution to (A.4.5) is: 
1 
Fe oes 2) 
ae Trae COTE 2 4 det C) (ee) 
Solvine for bi and c in (A.4.2) yields: 
ye 
Co Cle GF Cir 1C =a denne na 
ee eee Ee (A.4.7) 
1 
Ce ee PRCT ERC) 94 detrc)-, 
ee ee aa (A.4.8) 
Z 
Since I, = Tr C from (A.3.10): 
I_ =Tr cP (A.4.9a) 
ager ee (A.4.9b) 
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and the degree of polarization P can be written: 


I Pp 5 
depres hohe s 4 detae. 
P = Late C eo (A.4.10) 


Although the four real elements of an arbitrary coherency matrix 


Ke >) can be determined with several combinations of 


Cyy? xy!» 
polarizer and compensator readings, the following method is the most 


accepted. 


A.5 Experimental Determination of P | 
Rewriting equation (A.4.1) in terms of (A.3.12) and (A.3.13): 


-id 


I if 0 cos*x e siny cosx 
oi +1 (51) 


0 Ik ee adny cosx sin7x 


If a compensator is used to shift the phase of Hey relative to E,. by an 


amount 6 = -9, C’ remains the same and CP changes to: 
cos*x sinx cosx 
cP = I, (KES) 
sinx cosx sin2x 


This represents plane polarized light where ¢ = 0 or 7. C may then be 


written: 


I ai 0 cosy sinx cosx 
C Sigs +1 GAc53) 


0 M sinx cosx -sin*x 


By inserting a polarizer in the light with its pass direction oriented 


at an angle a with the x axis, the intensity can be described from 
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equation (A.3.5): 


It 
I(a) Sa + is cos“ (x - a) (A.5.4) 


When a = yxy, I(a) becomes a maximum: 
ps I (A.5.5) 


Conversely, when a = x + 7/2, I(a) is a minimum: 


il 


u 
Tin Teor (A.5.6) 


Finally, the degree of polarization can be written: 


M5 ea? tie 
P= 2 ee ee ae (Ax) 
qy Bs I) ae wi ain 


The only experimental difficulty using equation (A.5.7) is 
imposing a phase shift 6 = @ in the original light. As explained in 
Section A.2 this can be done with a quarter-wave plate. The major axis 
of the elliptically polarized part is first located by rotating a 
polarizer in the incident light. Then, by placing a quarter-wave plate 
along the direction of the major axis the elliptically polarized 
fraction is converted to plane polarized. The maximum and minimum 
intensities of the emerging light are determined with a polarizer and the 
degree of polarization is calculated from (A.5.7). In some cases the 
directions of maximum intensity are known a prtort. Equation (A.5.7) is 


then directly applicable without the need of a compensator. 
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APPENDIX B 


FLUORESCENCE 


B.1 Fluorescence 

As a consequence of the energy differences between electronic 
states certain molecules in solution will absorb electromagnetic 
radiation from the visible and ultraviolet parts of the spectrum. 
Several processes compete for the excess energy of an excited molecule 
and the relative kinetics of these processes are determined by the 
electronic structure of the molecule and its interaction with the 
surrounding medium. 

At room temperature most molecules occupy the lowest vibrational 
level of the ground electronic state and by absorbing light can be 
excited to higher vibrational levels of secondary electronic states. 
Sometimes the absorbed energy is sufficient to cause photochemical 
reactions such as dissociation but in most cases the excess enrgy is 
lost by emission or by collision with other molecules. 

Two electrons occupying the same ground state molecular orbital 
must satisfy the Pauli exclusion principle and have opposite spin. 
Hence, the total angular momentum from electron spin in this state is 
zero. In the excited state, however, an electron raised to an upper 
level can have the same spin alignment as the one remaining in the 
original orbital. This configuration is called the triplet state 
because there are three possible angular momenta associated with the 
total electron spin. Accordingly, the triplet state has a spin 


multiplicity of three and the ground state, or singlet, has a spin 


multiplicity of one. 
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There are two radiative processes by which an electronically 
excited molecule can return to the ground state: fluorescence and 
phosphorescence. 

Fluorescence is the most probable of these and is defined as the 
radiative transition between two states of the same multiplicity. The 
most common fluorescent transitions are from the first excited singlet 
to the ground state singlet and the average lifetime between absorption 
and fluorescence is approximately 1072 seconds. 

Phosphorescence is a more intricate process than fluorescence and 
is defined as the radiative transition between two states of different 
multiplicity. Such transitions are theoretically forbidden and only occur 
as the result of spin-orbit coupling (33). Because of spin-orbit coupling 
a molecule in a low vibrational level of an excited singlet can transfer 
to a high vibrational level of the corresponding excited triplet. This 
process is known as intersystem crossing (33). Spin-orbit coupling then 
enables the excited triplet to decay to the ground state singlet with the 
emission of light. Since the transition probability is very low, the 
lifetime of phosphorescence is rarely less than 10~* seconds and can be 
as much as one second. Phosphorescence is rarely observed from solutions 
because of various quenching mechant eas and can only be obtained from 
rigid media (49). The combined processes of fluorescence and 
phosphorescence are referred to as luminescence. 

During the lifetime of the excited state molecules in solution 
undergo many collisions and rapidly lose their excess vibrational energy. 
As a result the molecule relaxes to the lowest vibrational level of the 
excited singlet and the wavelength of subsequent emission is longer than 


that of the light originally absorbed. This feature is prominent in the 
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88 
emission spectra of solutions and is referred to as Stoke's law (33). 

Another important process competing for the energy of an excited 
molecule is internal conversion. This occurs when a molecule transfers 
from a low vibrational level of an excited state to a high vibrational 
level of a lower state of the same multiplicity. Again, the excess 
energy is lost in collisions and the process is non radiative. 

The main features of a luminescing molecule are shown in Figure 19. 
The diagram is drawn for diatomic molecules but is generally applicable. 
Due to spin-orbit coupling the triplet state of an atom consists of 
three closely spaced energy levels. These are not usually distinguished 
in molecules (33) and the triplet is illustrated as one level. In the 
triplet state there is less Coulomb repulsion between electrons and 
therefore this state lies below the corresponding singlet on a potential 
energy scale. 

From the Franck-Condon principle, light absorption takes place 
within the period of vibration of the light wave (= 10 1° sec) and during 
this time there is no appreciable nuclear motion. The transitions in the 
diagram are therefore drawn as vertical lines from the most probable 
nuclear configurations. In the scale of the diagram the energy spacing 
of the rotational levels is too smart to be shown clearly and only the 
vibrational levels have been drawn. These are represented with horizontal 
lines. 

The shape of the potential energy curve of an electronic level to 
which a transition is taking place determines the general features of an 
absorption or emission spectrum. The vertical lines in Figure 19 


represent the most probable transitions and correspond to the wavelengths 
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FIGURE 19. Franck-Condon principle applied to the luminescence 

of a diatomic molecule . The electronic levels shown are the ground 
¢ 

state singlet (S) , the first excited singlet (S ) , and the first 


4 
excited triplet (T ) 
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at which absorption or emission is a maximum. For absorption, the 
potential energy curve of the excited singlet (S”) is steeper to the left 
of the most probable transition and therefore the absorption spectrum will 
be broader to the shorter wavelength side of the maximum. If the 
vibrational energy spacing of the excited singlet is similar to that of 
the ground state singlet, the opposite is true for fluorescence and the 
emission spectrum is a mirror image of the absorption spectrum. 

Unless a luminescing system is continually excited the intensity 
of fluorescence decays exponentially. This follows from the fact that 
the number of molecules emitting per unit time is proportional to the 


number excited at that time; i.e., 


Son = ken (S211) 


where n is the number of excited molecules at any time. Letting Ny be 


the number of excited molecules at t = O and integrating yields 


n = nvee CB 2) 


Q cm Tare = kenp-e Celts 3a)) 
or defining k = 1/t: 
ner 
Diet ee (B.1.3b) 


where t is the mean radiative lifetime of fluorescence. 
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The above formulation assumes that every excited molecule under- 
goes fluorescence. In general, because of the many processes competing 
with fluorescence a fluorescence yield, das must be defined. 


by . Ginsteins of light emitted as fluorescence (B.1.4) 


einsteins of light absorbed 


where one einstein is an Avogadros number of photons. This does not 
affect the basic assumption of equation (B.1.1) and the decay of 


fluorescence remains exponential. 


B.2 Fluorescence and Electronic Structure 

Fluorescence is nearly always associated with the tm electron 
system of an unsaturated molecule (33). In the theory of atomic orbitals 
the electronic transition of lowest energy usually corresponds to the 
excitation of-an electron from a 7 orbital (bonding) to an upper 1* 
orbital (antibonding). However, for some molecules containing atoms with 
lone pair electrons, the transition of lowest energy corresponds to 
exciting an electron from a lone pair orbital to a m* orbital. These 
transitions have been classified by Kasha (27) according to whether the 
lone pair electrons conjugate fully, partially, or not at all with the 1 
electron system to which an Seeeeane being excited. If the lone pair 
orbital is located perpendicular to the nodal plane of the nm electron 
system then full conjugation is possible and the lowest energy transition 
is 1 > 1m*, For lone pair orbitals located parallel to the nodal plane no 
conjugation is possible and the lowest energy transition is from the non- 
bonding orbital to a m* orbital written as n> 7*. Located at inter- 


mediate angles, electrons from a lone pair can partially conjugate and 
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a7 
are only partially bonding. This absolves any strict planarity with 
respect to 7m electron configuration and the lowest energy transition is 
from the lone pair orbital to an antibonding orbital of nm origin written 
ae reg aie 

Fluorescence associated with these absorption transitions is 
strongest from molecules excited by 1 > m* and 1 > a. mechanisms. 
Molecules excited by n > m* transitions show relatively weak fluorescence. 
According to Wehry (49), there are two reasons for this. For (n,m1*) 
excited states the radiative lifetime is usually longer than for (1,17*) 
states and therefore nonavediarive processes competing with fluorescence 
such as internal conversion are more pip ieicanc: Also, since the energy 
differences between the first excited singlet and corresponding triplet 
are considerably less for (n,1*) states compared to (1,1*) states, 
intersystem crossing and phosphorescence are more probable. 

Because of their high degree of tm orbital conjugation most 
unsubstituted aromatic hydrocarbon solutions exhibit intense fluores- 
cence yield increases and fluorescence energies decrease as the length 
of a conjugated system increases. The most intensely fluorescent 
aromatic molecules are characterized by rigid, planar structures. The 
principal effect of increasing molecular rigidity is to decrease 
vibrational amplitudes which reduces the efficiency of intersystem 
crossing and internal conversion and thus increases the fluorescent 
yield (49). Figure 20 relates electronic structure and fluorescence for 


several aromatic molecules. 
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FIGURE 20. Fluorescence and electronic structure. 
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B.3 External Effects on Fluorescence 

When light absorption takes place, a molecule immediately 
(10719 sec) possesses an excited, geometrically extended electronic 
distribution. Because of the high inertia of the constituent nuclei the 
nuclear configuration cannot change in 10-!5 seconds and remains 
essentially the same as that of the ground level. This transitory 
structure is called the Franck-Condon excited state. During its 
radiative lifetime (10-9 sec) the excited molecule subsequently relaxes 
to an extended more stable nuclear geometry. Simultaneously the 
surrounding solvent molecules relax to a new configuration consistent 
with the electronic distribution and the geometry of the excited 
molecule. This "equilibrium" excited state is considerably lower in 
energy than the Franck-Condon excited state and all fluorescence 
originates from this level. Upon emission, the excited electronic 
distribution contracts to the ground state distribution leaving the 
nuclei in the excited configuration. This intermediate structure is 
called the Franck-Condon ground state. The molecule then undergoes 
structural and solvent reorientation and drops to the equilibrium ground 
state from which it originally absorbed. The energy relations of the 
four states are shown in Figure 21. 

External effects on the fluorescence of a molecule are determined 
by the chemical properties of the equilibrium excited state. 

Because the electron distribution in the excited state is less 
rigidly bound than in the ground state the excited molecule is usually 
more polarisable and therefore more easily solvated. Thus, by increasing 


solvent polarity the stability of (1,1*) excited states increases relative 


ve 


weasnibonas ohide Sour Fe Ry ois bef | olen | 

‘StHo boots besheriay a K [ Lento rt 

gia telouge tmags Seria) aid) ie eistagas tata ws | 

Ro hamet bas “bnoaae ra cs " anni 4 

riptesey art seas booets ana be 3) 
Ags we yut dase rh ese ire 

143, ldots “etapa, pa, 

tee Ness rit, een wae 6! ia: apfea 

ae aly 20 We ere Ts) uns aga 4 

uk wewo! Cile sobs 6a wt sande Renae’ ‘sont 

eons exobl t ae bre ode bysikoes: , re 

alyossugls Mwai es aid ne me 

sh cued suite Mathie tas Be | 

Ssoytabriu cael; vligaelihg oni) ao 

“Gowony aur disap 7] ite a ayo 1h. Sons di 


aris ae aitintios natn ott 


a “a 
| eee ‘fh a mia me  ennene 
r ant EAA YA, J ih | oa 


Pontes 96 amb ah aacrsviled to ieee wie: Ps 
‘Otain. ja sada suv ope aly Ye. | 

‘gst et oe a s Pi Bit ch Be dudeiankh é ° 
ef Lise Ri a: 1-39 Loin by sh koe” a $a 10 | ponte 


. re * okt ae Hh 
é ih 


sonore 8. and -be1sitton aoe 


Vani sey é ; -* 


a Ti 


FIGURE 21 . 


95 


FRANCK-CONDON 
EXCITED STATE 


EQUILIBRIUM 
EXCITED STATE 


ABSORPTION 
FLUORESCENCE 


FRANCK-CONDON 
GROUND STATE 


EQUILIBRIUM 
GROUND STATE 


Effect of solvation on the absorption and fluorescence 


of a luminescent molecule. 
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to the ground state and m* + 7 fluorescence shifts to longer wavelengths. 

For molecules containing highly electronegative atoms such as 
oxygen or nitrogen, the ground state is more easily solvated. In the 
excited state the charge density around the electronegative atoms is 
lowered and the resulting dipole moment is considerably less than in the 
ground state. Increasing the hydrogen bond capacity or the polarity of 
the solvent in which these molecules are dissolved causes the stability 
of the ground state to increase relative to the excited state and there- 
fore the absorption wavelength decreases. Since (n,7*) lowest excited 
states are commonly associated with fluorescent molecules containing 
nitrogen or oxygen, n > 7T* absorption incurs a blue shift as solvent 
polarity is increased. 

The influence of external pH on the fluorescence of a molecule 
containing acidic or basic functional groups is often drastically 
different from that observed for its absorption spectrum. In order for 
acid-base reactions in the excited state to perturb fluorescence spectra 
it is necessary for protolytic dissociation to be appreciably more rapid 
than the rate constants for fluorescence. Brédnsted acidity differences 
between ground and excited states of organic molecules are large, 


ranging from 4 to 9 pK, units (49).! As an example of the effect of 


IThe pH of a solution is defined as the negative logarithm (base 
10) of the hydrogen ion concentration; i.e., pH = - log[H’]. For one 
liter of water at 23° C the dissociation EXO — di + OH- produces a Ht 
concentration of 10-’ M/1l. Thus, neutral pH is 7 and pH values less than 
or greater than 7 are acidic or basic respectively. 


According to the law of mass action the dissociation of a weak 
acid can be described by an equilibrium constant; 
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external pH, fluorescein exhibits a brilliant yellow-green fluorescence 
in alkaline solutions and a weaker blue-green fluorescence in acidic 
are (36)).. 

There are two important types of processes by which a solvent or 
extraneous solute can quench the fluorescence of a molecule in solution. 

The first type of process is an interaction between solute and 
solvent that increases the rates of radiationless processes in the 
excited state. Most noteworthy of these is the quenching of fluorescence 
by molecular oxygen. Because 05 is paramagnetic in the ground state 
(i.e., triplet ground state) it can induce intersystem crossing in 
excited solutes and thus decreases the fluorescence yield. 

The second of these processes is the direct transfer of the 
electronic energy from solute to quencher. Direct transfer occurs in 
one of three ways. If the absorption spectrum of the solvent overlaps 
in the fluorescence spectrum of the solute then it is possible for the 
solvent to quench fluorescence by direct absorption. However, unless 
conditions are ideal this will not occur to any great extent. A more 
important process involving direct energy transfer is collisional 
transfer. Collisions are diffusion-controlled and the efficiency of 
collisional transfer depends on Ene iseosicy of the solvent. The third 
direct transfer mechanism is known as "long range resonant transfer" or 
Forster transfer (33). This is a non-radiative intermolecular dipole- 


dipole transition that does not require close approach of donor and 


where iA te Ag. Lhe pK,. of an acid is the negative logarithm (base 
10) of Kg and thus represents the degree of dissociation. 
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acceptor. The efficiency of this mechanism can be expressed in terms of 
a critical separation distance at which the probability of transfer is 
equal to that of spontaneous decay. This process is .only efficient when 
there is large overlap between emission and absorption spectra of donor 
and acceptor molecules respectively, and when the fluorescence yield of 
the donor is high. 

The last external influence treated in this section is temperature. 
In the absence of quenchers, fluorescence efficiency is usually a maximum 
around room temperature and decreases slowly with increasing temperature. 
This is because the probability of internal conversion and intersystem 
crossing increases at higher temperatures which reduces the fluorescence 
yield. At low temperature the probability of these processes reaches a 
limiting value corresponding to internal conversion and intersystem 
crossing from the lowest vibrational level of 3“. The variation of 


fluorescence yield can frequently be represented by the equation (4): 


1 e7E/RT 
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APPENDIX C 


POLARIZED FLUORESCENCE 


C.1 Polarized Fluorescence 

As discussed in Appendix B.2 absorption and fluorescence in the 
visible or near ultraviolet are usually associated with m electron 
systems. Electron conjugation in most molecules defines fixed axes along 
which absorption and fluorescence preferentially occur and such molecules 
may be treated as electric dipole oscillators (19). 

The most general classical model of a fluorescent molecule 
consists of an emission oscillator E and an absorption oscillator A which 
have different frequencies but are coupled together so that light absorbed 
along A is reemitted along E. The directions of E and A can coincide but 
generally intersect at some angle 8. Polarized fluorescence is a direct 
consequence of totally anisotropic (linear) oscillators. 

If the emission axes from a group of excited molecules are randomly 
oriented in solution it follows from symmetry that the subsequent 
fluorescence will be unpolarized. Thus, in order to produce polarized 
fluorescence it is necessary that the excited molecules retain a non- 
random orientation relative to a fixed axis. This is normally a difficult 
requirement because during the lifetime of the excited state molecular 
rotation and energy transfer tend to depolarize fluorescence. However, 
by using a sufficiently viscous solvent and a low enough solute 
concentration these effects can be minimized (Sections C.2 and C.3). 

When depolarization from rotation and energy transfer are negligible the 


degree of polarization of fluorescence depends on three factors; the 
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method of excitation, the direction of observation, and the coupling 
angle between A and E in the individual molecules. 

The simplest system to analyze is shown in Figure 22 where it is 
assumed that depolarizing effects are unimportant and where A and E 
coincide. If this system is irradiated with vertically polarized light 
(Ey along OZ) then for any oscillator the projection of A along Eo is 
proportional to cos@ and the probability of excitation is equal to cos“6. 
Absorption preferentially occurs when A is parallel or nearly parallel to 
the Z axis and since A and E coincide the same cos70 dependence applies 
to fluorescence. Thus for a given number of excited molecules the 
fluorescence intensity viewed through a plane polarizing filter will be 
a maximum when the pass direction is oriented parallel to the Z axis and 
a minimum when oriented perpendicular to the Z axis. This is shown in 
Figure 22 where the Y axis has been chosen as the direction of observa- 
tion. For vertically polarized excitation the Z axis is an axis of 
symmetry and viewed along any direction in the XY plane I, and I, remain 
constant. Equation (A.5.7) is directly applicable and the degree of 


polarization of fluorescence is given by: 


Laxy es nine eee et 
ee ere RSC ee EE plas 
FT Sisco TSE TREE Neigh ct) 
max min 
where the subscripts n and 4s refer to directions parallel and 


perpendicular to the axis of symmetry defined by Eo of the exciting light. 
When equation (A.3.10) is extended to three dimensions the total emitted 
intensity can be described as the sum of three orthogonal components and 


because of the symmetry resulting from vertically polarized excitation 
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VERTICALLY POLARIZED 
EXCITATION 


FIGURE 22 . Absorption and emission of a randomly oriented electric 


dipole oscillator. Absorption and emission axes coincide along OE . If 
the direction of observation is along the Y axis , the fluorescent 
intensity viewed through a plane polarizing filter will vary with the 


orientation of the pass direction @ as shown. 
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[eee Te Or | (C212) 


where Iy is the component directed along the Z axis and the two equal 
I, components are in the XY plane. 

Equation (C.1.1) employs only one of the I, components because the 
remaining component is parallel to the direction of observation. 
However, if the direction of observation is inclined at an angle y to the 
XY plane (Figure 23) ae becomes the sum of the projections of Iy and 
the second I, component along a direction perpendicular to OP and in the 
plane OPY”. As y increases the projection of Iy is reduced by the 
factor cosy and that of I, is increased by the factor sin’y. Me can 


therefore be written: 
I (y) = Iy cos?y + 1 sin?y (Cal 3) 
I is also perpendicular to OP but lies in the XY plane and is there- 


min 


fore independent of y. From symmetry: 
Line) = yl CG te24)) 
and the degree of polarization of fluorescence becomes: 


(Iy cos*y + I, sin*y) - ly 
XG) Gh ae ee ee CGAL) 
(Iy cosy te I sin?y) ce yy 
From equation (C.1.5) P is a maximum when viewed in the XY plane (y = 0) 
and decreases with increasing y until the direction of observation 


coincides with the Z axis and P = 0. Thus for vertically polarized 


excitation P is normally measured along the Y axis as in Figure 22. 
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VERTICALLY POLARIZED 


EXCITATION 

E 

0 
(a) 

UNPOLARIZED EXCITATION 
Ty 
(b) 
Ty 
FIGURE 23. Fluorescent intensity components and symmetries from 


systems excited by (a) vertically polarized light directed along 


the X axis (b) unpolarized light directed along the X axis 
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Polarized fluorescence can also be produced by excitation with 
unpolarized light. If the same system is excited by a beam of unpolarized 
light directed along OX then with respect to absorption and fluorescence 
the X axis will become an axis of symmetry. Equation (C.1.2) can again 
be applied but for unpolarized excitation the two I, components lie in 
the ZY plane while the Iy component lies along the X axis (Figure 23). 
Because there are two unequal components in the XY plane the degree of 
polarization does not follow the same format of equation (C.1.5). The 
method for calculating P is the same however and viewed along OP as in 
Figure 23: 
[ (I, cos? = sin?¢) sin?y + I, cos*y] - (Ty cos2@ + I, sin?) 
aye [(I,cos* + I, sin?) sin?y + I, cos*y] + G, cos? + I, sin?) 
(G5L26) 
where the subscripts refer to directions relative to the axis of symmetry. 
When ¢ = y = O the direction of observation coincides with the Y axis and 
equation (C.1.6) reduces to: 


Tysseey 


a =e GL. 
I, + Iy ( ) 


This is the maximum value of P obtainable from equation (C.1.6) and 
therefore for unpolarized excitation, the Y axis is the preferred 
direction of observation. Although unpolarized excitation is of 
theoretical interest the degree of polarization of fjuorescence is always 
greater from systems excited with plane polarized light and for this 
reason plane polarized excitation is normally used. 

In the absence of depolarizing effects and for a given method of 


excitation, the maximum polarization that may be produced from a system 
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is referred to as the principal polarization and is designated by Po: 
Referring again to Figure 22, the vertical and horizontal intensity 


components from a single oscillator are given by: 


Iy & cos“6 (Ce 8) 


I, «< sin*6-sin76 (GC. 5.9) 


If equations (C.1.8) and (C.1.9) are weighted to include the probability 
of excitation and integrated over 6 and 9, Po can be determined from 
equations (C.1.1) and (C.1.7). For systems irradiated with vertically 
polarized light the probability of excitation varies as cos*@ and 

a nae I,)/(1y + 1.) = 1/2. If the exciting light in unpolarized 
PO 
be obtained by a less formal method. Weber (47) has shown that for n 


(ip =i, )/ CG; +i, ) = 1/3.) Theses limiting polarizations can also 


individual oscillators, P may be expressed by: 
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where the plus sign refers to unpolarized excitation and the minus sign 
refers to vertically polarized excitation. In equation (C.1.10) Ps is 
the polarization from the "ith oscillator and f£, is the contribution from 


the "i'th oscillator to the total fluorescent intensity, i.e.: 


n 
f; = I,/l.; where Rey £5 = 61 Geleid) 


Since f can represent the probability of excitation and P can be 


calculated for a given 9 and $, equation (C.1.10) can be summed over over 
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all 6 and ¢ to give 1/2 or 1/3 depending on the method of excitation. 
The principal polarizations of 1/2 and 1/3 apply only to systems 
containing molecules whose absorption and emission axes coincide. 
Jablonski (24) has shown that if the angle between the absorption and 


emission oscillators is 8, the principal polarization can be described 


by: 


ne cos*B - 1)/(cos28 + 3) (Cade 2) 


for vertically polarized excitation, and 


Po = (3 cos2g =4)/7 - cos2f) CCR ee 3,) 


for unpolarized excitation. Because 8 can only vary between O and 71/2, 
P will lie between the limits of +1/2 and -1/3 for equation (C.1.12) or 
+1/3 and -1/7 for equation (CEng. 

The importance of equations (C.1.12) and (C.1.13) is experimentally 
demonstrated by the variation of the degree of polarization of 


fluorescence with exciting wavelength. Figure 24 shows the absorption 


and polarization spectra of a basic glycerol solution of fluorescein. ! 


lThe polarization spectrum, or polarized excitation spectrum, 
shown in Figure 24 was constructed with data obtained from a Perkin-Elmer 
MPF-4 spectrophotometer. A description of the instrument is given in 
Chapter Two. To obtain a polarized excitation spectrum, a sample is 
irradiated with plane polarized light and the fluorescence intensities 
parallel and perpendicular to the electric vector of the exciting light 
are analyzed through a plane polarizing filter. The fluorescence 
components are observed at a fixed emission wavelength, em, while the 
exciting wavelength is varied from the near ultraviolet through to the 
visible. The polarization at any exciting wavelength is given by: 
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The various peaks in the polarization spectrum correspond to different 
absorption transitions within the molecule and in accordance with 
vertically polarized excitation the curve varies between the limits of 
+1/2 and -1/3. These limiting polarizations cannot be reached 
experimentally however, and as shown in Figure 24 the maximum and minimum 
polarizations are .43 and -.23 respectively. Feofilov (18) contends that 
this is due to incomplete anisotropy in the oscillators and represents a 
breakdown of the classical model in completely describing polarization 
phenomenon. Because of incomplete Oe eeee Po can never be accurately 
determined and therefore equation (C.1.12) has limited experimental 
application. Feofilov has alleviated this problem somewhat by deriving 
an expression which relates various limiting polarizations to molecular 
symmetry (18). With this formulation many polarization spectra can be 


resolved. 


C.2 Depolarization by Energy Transfer 


In sufficiently viscous solutions where rotational depolarization 
is minimal Feofilov and Sveshnikoff (19) empirically found that P varies 
inversely with solute concentration, i.e.: 

1 
= —— + AeCeT (Cy Ze) 


PO 


rd|h 
| 


where t is the lifetime of the excited state, c is the solute 


emission monochromator. 


For the fluorescein solution, depolarization effects were minimal, 
and Figure 24 depicts the variation of the principal polarization (P,) 
with the exciting wavelength. The wavelength of maximum emission was 
Nem = 515 nm, and rex was varied from 250 nm to 550 nm. At an emission 
bandwidth of 10 nm (slit width = 3 mm), G was calculated from the method 
of Azumi and Mclynn (1) to be G = .454. This spectrum agrees with that 
previously published by Chen and Bowman (17). 
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FIGURE 24 . Absorption and polarization spectra of fluorescein , 


Dex 10° M in 95% glycerol and 5% 0.10 N NaOH by volume, recorded 
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concentration and Po is the principal polarization. 

As an example, the polarization from basic glycerol solutions of 
fluorescein begins to decrease at concentrations > .3 mM and falls to 0 
at approximately .03 M. At these concentrations the average inter- 
molecular distance between solute molecules is about three times their 
kinetic radii and energy transfer can only be explained by Férster 
transfer (18, Appendix B.3). By considering the probability of long- 


range energy transfer between two dipole oscillators Weber (45) has 


derived a concentration-dependent polarization. According to Weber: 


eRe Oe -3 

1 : 1 bak NN , 4teN vA 10 | 

roeting pe art ees 
) 15° (2a) 


Pee (Ce2 2) 


where the plus sign refers to unpolarized excitation and the minus sign 
refers to vertically polarized excitation. In equation (C.2.2) Py is 
the principal polarization, N is Avogadro's number, Ro is the critical 
radius between parallel dipoles at which the probability of emission 
equals the probability of transfer, 2a is the molecular diameter, and c 
is the solute concentration in moles/liter. 

Applying this result to the data of Feofilov and Sveshnikoff, 
Weber has calculated a critical radius of 26 A for fluorescein. By 
keeping the concentration sufficiently low that the average distance 
between solute molecules is greater than the critical radius, concentra- 


tion-dependent polarization can be minimized. 


Gese Rotatzonal! Depolarization 


When concentration effects are negligible depolarization of 


fluorescence depends only on the interval between excitation and emission 
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and on the speed of molecular rotation. Various models have been 
formulated to describe rotational depolarization and all have been 
derived from the classical theory of Brownian motion. This theory assumes 
that equipartition of energy (kT/2 average rotational energy per degree of 
ee eed ot) applies to the motion of the molecules and that molecular inter- 
action can be described by macroscopic properties. The first analytical 
description of rotational depolarization was done for spherical molecules 
by Perrin (35). This was later revised by Soleillet (42), and more 
recently Weber (48) has extended the theory to include ellipsoidal 
mollecular geometries. 

Following the description by Weber, the angular speed of a molecule 


rotating in a viscous medium is given by: 


W = Wo e (Fr/I)t C62 321) 


where Wo is the initial angular speed, I is the moment of inertia of the 
molecule, and 1 is the damping coefficient which is a function of n the 
viscosity of the medium and the geometry of the molecule. From equation 
(C.3.1) it can be shown that the average angle rotated by the molecule in 


atime intenva let is): 
<sN0e >o= he < Bo At / fe (C32) 


where < > indicates time average and where <E> is the average rotational 
energy of the molecule. Equation (C.3.2) is a fundamental property of 
Brownian motion. 

The diagram below shows two successive rotations of a molecule 


about independent axes. The position of a direction rigidly bound to a 
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sphere lies successively along OA, OB and OC making angles a, 8 and y 


as shown. 


y is the azimuth of OC about the plane AOB 


From elementary geometry: 
cos Y = cos a-cos 8 + sin arsin B-cos p (Cy 32 3) 


Averaging equation (C.3.3) for a large number of molecules that have 


undergone two successive rotations: 


<cos Paes = <cos a> ,*<cos B>. ap <isilin a> +<sin B>,*<cos p> 


B v 


(Cz34) 
If the molecules undergo isotropic rotation all values of y are equally 
probable so that <cos Pa: = 0 and <cos*p> = 1/2. Thus, by squaring 


equation (C.3.4): 
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2 = 2 5 2 Paks ROS ad Se LSS 
<cos Lie <cos Ones <cos B>, ae ; (C255) 
or by trigonometric rearrangement: 
3 <cos*y>_ - 1 3} <cos*a> - 1 3 <cos*B>, - 1 
ae ee ee er eee gh sete a et aes (G36) 
Z 2 2 


Equation (C.3.6) describes y as the product of two angular impulses and 
the effect of any number of individual impulses may be computed in the 


same way. 


For a small rotation a, cosa = 1 - sinta ~ 1 - a%, and using 


equation: (G. 3.2) 


costa = 1 - 4e<E>-At/f£ (CA3e7) 
Dawe 
Se Asosia at = 1 - 6+<E>-At/£, (C.3.8) 


The value of <cos*y>_ at any time t may be calculated by dividing t into 


an arbitrary number of small intervals and applying equations (C.3.6) and 


(C7338). ~Fhus, 


3 <cos*y(t)>, - 1 
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Equation (C.3.10) was derived for isotropic rotation and in general is 
only applicable to spherical molecules. Since there are two degrees of 
freedom for the rotational motion of a sphere about two independent axes, 
the average rotational energy is given by <E> = kT, where k is Boltzmann's 
constant and T is the absolute temperature. In addition, the damping 
coefficient for a sphere rotating in a viscous medium is given from 
basic hydrodynamics (Kirchoff) as f. = 8emener, where r is the radius of 
the sphere and n is the viscosity of the medium. Substituting these 


values in equation (C.3.10) and defining: 
Py = 4emener3/k-T (ager) 
as the rotational relaxation time, 
“ees ( 0) = 1/3 + 2/3 eo 3t/P0 (G3 502) 


The concept of isotropic rotation is useful because if the average angle 
of the emission oscillator with a given axis is known at t = 0 the average 
angle at any further time can be calculated from equation (C.3.6). 

From Figure 22 the fluorescent intensity components from a single 
oscillator were given by equations (C.1.8) and (C.1.9) as Iy = cos*6 
and I,=sin76+sin*9. For systems irradiated with vertically polarized 
light all values of $ are equally probable and <sin*>, = 1/2. TAUSs, Lor 
all oscillators making an angle 9 with the Z axis the contribution to the 


total perpendicular component is given by: 


a8 
iy ee cea (¢.3.13) 


and the polarization from these oscillators becomes: 
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Ty q, 2. 
30C0sc0) =—1 i 
Ip CP SS iS Vere EE ees GC25- 4) 
ae ey eae cos?6 +1 
or 
Li pee Pe 2/3 


Z 


Using equation (C.1.10), the polarization from all oscillators making 


variable angles 6 is: 


1 1 a! 1 
SS ee (C3550) 
y : Gan /2 £(6) O=n/2 3 cos*6 - 1 
z z £(6)° 
1 Al 2 
6=0 poy 8=0 
e] 
6=1/2 
and since 2& £Copvrasr: 
6=0 
LP Ae ommbedt et 2/ Bier (Casa) 
ie 3 3 <cos20ea = 1 ee 
22 


The temporal dependence of equation (C.3.17) can be determined by 


considering the separate polarizations at times 0 and t after excitation. 


At t = 0: 
Leelee Pall (G23 718) 
Po Dur 3 <cos*6(0)>, -1 
2 


and at time t after excitation: 


1 it 2/3 
eee CCr3ek9) 
Pie fo ce <cos*6(t)>g - 1 
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If the rotations are isotropic: 


3 <cos*6(t)>g = legs <cos*6(0)>¢ Seek gt We S <cos*a(t)>, - 1 
eae Ne Oe ee a a 
(Og o20) 
where a(t) is the angle swept by the emission axes between times O and t, 
and by substituting equation (C.3.20) into (C.3.19): 


pe Se ee Gene APNG 2) (ae Geka ae We 2 eee Bae (Gear 2) 
P ) 3 <cos*a(t) y> - 1 


Although equation (C.3.21) expresses P as a function of time, it 
does not take into account that the number of molecules emitting at any 
time are distributed according to when they were excited and thus 
according to the average angle they have rotated relative to their 
original positions. To calculate the polarization observed at any time 


after excitation, equation (C.1.10) can be used again giving: 


eR ia eg et Sel eee le 
P 3. «t= P 3 t=00 
y EGE 0 ee <cos*a(t)> *f(t) -1 
1 1 3 fe} 
ED Sn 5 t=0 
Pe 
(Ca3n22) 


where f(t) is the fraction of the total intensity due to molecules 
emitting t seconds after excitation. The function f(t) is given by 


equation (B.1.3) and represents a constant rate of emission, i.e., 


entit (O83 223) 
EF 


Pinte 


where t is the lifetime of the excited state. In the denominator of 
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equation "(C.3% 22): 


t=~ co es 
De <cos“a(t)> 9° f(t) as id tos-a(b)> = 
t=0 0 


dt (Ci5624) 


and since <cos*a(t)> = 1/3 + 2/3 e 3t/Pg for spherical molecules, 
equation (C.3.24) can be integrated and inserted in equation (C.3.22) 


giving: 


1 iL ENG 
= SS SES ° ale + — Gs aD 
| ) Po (CES25)) 


Thus, derived from the classical theory of Brownian motion, rotational 
depolarization depends on the ratio of the average excited lifetime and 
the rotational relaxation time of the molecular oscillators. 

Equation (C.3.25) has been derived for vertically polarized 


excitation and a similar derivation for unpolarized excitation yields: 


teeta tla es | ellen, 2s tw Lap uric (C.3.26) 
Bie 3 dy 


As they appear, equations (C.3.25) and (C.3.26) are difficult to verify 
experimentally. Since the volume of a molecular sphere of radius r is 


given by V = 4/3nx°, equation (C.3.25) can be modified giving: 


td|h 


-< (3574): 1+, (G327) 


Perrin (35) and Pringsheim (36) have verified equation (C.3.27) for 
several solutions by obtaining a straight line for graphs of 1/P against 


T/n. Although this method is somewhat limited in that it assumes V and 
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t are independent of T/n, the experimental evidence corroborates the 


classical derivation of rotational depolarization. 
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APPENDIX D 


COMPUTER PROGRAM FOR DETERMINING Ee 


The computer program for calculating P. was written in Fortran 
and was executed by an Artronix PC-12 interactive system. The system 
compiler did not accept logical IF statements and therefore the 
arithmetic IF statement was used. PGCTL is a "system" subprogram for 
controlling the oscilloscope display and can be deleted for operation on 
other computers. A minimum memory of 16 K is required to run the 
program. 

The main features of the program are shown by a flow diagram in 
Figure 10. To provide a least squares fit to I, and I, data, a sub- 
routine CURF was used which is a modification of a program originally 
designed by Bevington (3). The important variables in the program are 


given in the following list. 


X = reaction time 


Y = Ty and Ty data 


G = grating correction factor 

S = standard deviation in the Iy, and I, points 

R = corrected ratio of and 14 background components 
BB = background component (Figure 8) 

TF = fictitious time of filtration 


C = least squares fit to quadratic polynomial for Iy data 
B = least squares fit to quadratic polynomial for I, data 


FB = starting filtrate (Te) for the search for optimum Ey 
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EL. 
DFB = AI 
MNM = number of iterations 
CLCD1 = subprogram which calculates Kp from a value of Le or 


e : ; 
I, at a particular time 


DB = Ky 

DC = Ky 

PT =Stotal’polarization P(t) 

PC = polarization from the cytoplasm P (Kt) 
PF = background polarization P(t) 


YET1 = equation (2.5.13) for Dace) and Loce) 


IC = R, from equation (2.4.4) 


PVM(1) = Ky 
PVM(2) = SEM? 
PVM(3) = 0 
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DIMENSION X(30) ,Y(30) ,SGMY (30) ,A(3) ,DLTA(3) ,SGMA(3) , YFT(30) 
DIMENSION X1(30),Y1(30) ,B(3),C(3) ,V1(8) ,V2(8) ,V3(8) ,V4(8) 


DIMENSION V5(8),PVM(15, 3) 


CALL PGCTL(1) 

READ(1,1)G,BB,R,S 

FORMAT (G"G 'G="\F10.4,"" BB=",F10.4," R=", F10.4," S="', F1O.4) 
READ (1,100) TF 

FORMAT(" "," TF=""\F10.4) 

READ (1,2)NPTS,NTRMS, MODE ,N 

FORMAT(" "," NPTS="",14," NTRMS="",14," MODE=",14," N=",14) 
BC = BB*R 

DO 114 I=1,2 

DO 3 IA=1,NPTS 

READ (1,4) X(IA) ,Y (IA) 

FORMAT(" ",") X="",E13.5," Y=",E13.5) 
IF(I-2)6,5,5 . 
Y(IA)=Y(IA)*G 

SGMY (IA) =S 

CONTINUE 

K=NPTS-N 

DO 7 IB=1,K 

X1 (IB) =X (IB4+N) 

Y1(IB)=Y(IB+N) 

CONTINUE 

DO 8 IC=1,NTRMS 

DLTA(IC)=.005 

CONTINUE 

A(1)=BB 

TF G2) 9210.10 

A(1)=BC 

A(2)=1.5 

A(3)=-.005 

FL=.001 

CALL CURF(X1,Y1,SGMY,K,NTRMS,MODE,A,DLTA, SGMA, FL, YFT, CHI) 
DO 12 ID=1,NTRMS 

WRITE(1,13)A(ID) 

FORMAT(" "," A="",E13.5) 

CONTINUE 

WRITE (1,14) CHI 

FORMAT(" ","" CHI=",E13.5) 

READ (1,15) IE 

FORMAT(" "," REFIT WITH NO. .GE. 1 ",14) 
IF(IE-1)16,11,11 

V=YTT1(A, TF) 

IX=NPTS+1 

X(IX)=TF 

Y (IX) =V 

FL=.001 , 

CALL CURF(X,Y,SGMY, IX,NTRMS,MODE,A,DLTA, SGMA, FL, YFT, CHI) 
WRITE(1, 21) CHI 

FORMAT(" "," CHI="",E13.5) 

DO 18 IG=1,NTRMS 
WRITE(1,19)A(IG) , SGMA(IG) 

FORMAT(" "," A="",E13.5," S=",E13.5) 
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CONTINUE 

READ (1,20) IH 

FORMAL CGH UIREFIT WITH NO...) )GE2 9 1) © "4 14) 
IF (IH-1).22,17,17 

DO"23) Li=1, 1X 

LP(I=2)24,25,25 

YFT (1I1)=YFT(1I1L)/G 


WRITE (1,26) YFT (IIL) 
FORMAT (S"*) 6" RYETH""$E13\55) 
CONTINUE 


DO 114 IJ=1,NTRMS 
IF (1-2) 28,29, 29 
C(IJ)=A(IJ) 

60 To 414 
B(IJ)=A(IJ) 
CONTINUE 


CALCULATION OF POLARIZED FLUORESCENCE 


READ (1,80)FB,DFB,MNM 
FORMAT(" "," FB=s",F10.4," DFB=",F10.4," MNM="",14) 
DO 150 IS=1,MNM 

PVM(IS,1)=FB 

FC=FB*R 

CALL CLCD1(B,FB,TF,DB) 

GALL CLCD ECC, FC.TF DC) 
SUMP=0. 

T=10. 

DO 81 IK=1,8 

PrsyTT.(b,1) 

P2=YET1(B,DB,T) 

P3=YTT1(C,T) 

P4=YETI(G.DC,.L) 
PT=(P1-P3) / (P1+P3) 
PC=((P1-P2)-(P3-P4) )/ ((P1-P2)+(P3-P4) ) 
PF=(P2-P4) /(P2+P4) , 
V1 (IK) =PT 

V2 (IK) =PC 

V3 (IK) =PF 

V4(IK)=1. - ((P2+P4)/(P1+P3)) 
SUMP=SUMP+V2 (IK) 

T=T+2. 

CONTINUE 

PA=SUMP/8. 

PVM(IS,2)=PA 

SUMS=0. 

DO 82 J=1,8 
V5(J)=ABS (PA-V2 (J) ) **2 | 
SUMS=SUMS+V5 (J) 

CONTINUE 

SP=SQRT (SUMS) 

PVM(IS,3)=SP 

FB=FB+DFB 

CONTINUE 
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WRITE(1,83) 

HORMAT Cua iO XS ERY LS Xs PAM 15%, SP") 

DO 151 IT=1,MNM 
WRITE(1,152)PVM(IT,1),PVM(IT,2),PVM(IT, 3) 

BORMAT Ctl) 5X, F1004 5X, F004, 5X, E13: 5) 

CONTINUE 

PAUSE 

READ (1,84) IL, IM 

PORMAT Ce! .') READLEDT: WETHEL 1490) VIEWSPTS WITH 1) >", 14) 
IF (IM-1) 85,88, 88 

T2110; 

WRITE(1,91) 

ROMANCE LT KOT CMLL Xe! Pie LIke PC rT Xe PR LT XetLE') 
DO 93 IN=1,8 

WRITE (1,90)T,V1(IN) ,V2(IN) ,V3(IN) ,V4 (IN) 

FORMAT Gui e2XRE TO, GE OX  FlOSA- 2k PLOnG 2k sF 0.4. 2X, 610.4) 
T=T+2. 

CONTINUE 

IF (IL-1) 89, 86, 86 

GO TO 99 

STOP 

END 


CURVE FITTING SUBROUTINE 


SUBROUTINE CURF(X,Y,SGMY,NPTS,NTRMS,MODE,A, DLTA,SGMA, FL, YFT,CHI) 
DIMENSION X(25),Y(25),SGMY (25) ,A(3) ,DLTA(3) ,SGMA(3) , YFT(25) 
DIMENSION WT (100) ,ALP(10,10) ,BET(10) ,DRIV(10) ,ARRAY(10,10) ,B(10) 
NFREE=NPTS-NTRMS 

IF (NFREE) 13,13, 20 

CHI=0. 

GOmTORLLO 

DO 30 I=1,NPTS 

IF(MODE)22,27,29 

HHGVGE)) Co ean 25 

WT (1)=1./Y(1) 

GO “£0: 30 

WT (1)=1./(-Y(1)) 

GO TO 30 

WI(I)=1. 

GO TO 30 

WT (1)=1. /SGMY(1)**2 

CONTINUE 

DO 34 J=1,NTRMS 

BET (J)=0. 

DO 34 K=1,J 

ALP (J ,K)=0. 

DO 50 I=1,NPTS 

CALL FDRIV(X,1,A,DLTA,NTRMS, DRIV) 

DO 46 J=1,NTRMS 

BET (J) =BET (J )+WT (1) * (Y (1) =QPOL (X,1,A)) *DRIV (J) 

DO 46 K=1,J 

ALP (J ,K)=ALP (J, K)+WT (1) *DRIV (J) *DRIV (K) 

CONTINUE 
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51 DO 53 J=1,NTRMS 
DO 53 K=1,J 
53 ALP (K,J)=ALP(J,K) 
61 DO 62 I=1,NPTS 
62 YFT (1)=QPOL(X,1,A) 
63 CHI1=FCHI (Y, SGMY ,NPTS ,NFREE, MODE, YFT) 
yak DO 74 J=1,NTRMS 
DO 73 K=1,NTRMS 
73 ARRAY (J, K)=ALP (J ,K) /SQRT (ALP (J, J) *ALP(K,K) ) 
74 ARRAY (J,J)=1.+FL 
80 CALL MINV(ARRAY, NTRMS, DET) 
81 DO 84 J=1,NTRMS 
B(J)=A(J) 
DO 84 K=1,NTRMS 
84 B(J)=B(J)+BET (K) *ARRAY (J ,K) /SQRT (ALP (J , J) *ALP (K,K)) 
91 DO 92 I=1,NPTS ; 
92 YFT(1I)=QPOL(X,1,B) 
93 CHI=FCHI (Y, SGMY ,NPTS , NFREE, MODE, YFT) 
LF (CHI1-CHI)95,101,101 
95 FL=10. *FL 


GOTO 71 
101 DO 103 J=1,NTRMS 
A(J)=B(J) 
103 SGMA(J)=SQRT(ARRAY(J,J)/ALP(J,J)) 
FL=FL/10. 
110 RETURN 
END 
EG 
C REDUCED CHI SQUARE FIT FUNCTION 
C 


FUNCTION FCHI(Y,SGMY,NPTS,NFREE, MODE, YFT) 
DIMENSION Y(25),SGMY(25) ,YFT(25) 

11 CHI=0. 

ity LF (NFREE) 13, 13,20 

13 FCHI=0. 
GO TO 40 

20 DO 30 I=1,NPTS 

at IF (MODE) 22,27,29 

22 LF(Y(1))25,27,23 

23 WT=1./Y(I) 


GORTO L380 
25 WT=1./(-Y(1)) 
GO TO 30 
27 WI=1. 
GO TO 30 


29 WT=1. /SGMY (1) **2 
30 CHI=CHI+WT* (Y(1)-YFT (I) )**2 
oH, FREE=NFREE 
32 FCHI=CHI/ FREE 
40 RETURN 
END 
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MATRIX INVERSION AND DETERMINANT SUBROUTINE 


SUBROUTINE MINV(ARRAY,NORD, DET) 
DIMENSION ARRAY(10,10),1IK(10),JK(10) 
DET=1. 

DO 100 K=1,NORD 

AMAX=0. 

DO 30 I=K,NORD 

DO 30 J=K,NORD 

IF (ABS (AMAX) -ABS (ARRAY (1,J)))24, 24, 30 
AMAX=ARRAY (I,J) 

IK(K)=1 

JKC) =I 

CONTINUE 

IF (AMAX)41,32,41 

DET=0. 

GO TO 140 

I=IF (K) 

PE@l=K)21,51,43 

DO 50 J=1,NORD 
SAVE=ARRAY (K, J) 

ARRAY (K, J) =ARRAY (I,J) 
ARRAY (1, J)=—-SAVE 

J=IK (K) 

IF (J-K) 21, 61,53 

DO 60 I=1,NORD 

SAVE=ARRAY (1,K) 

ARRAY (L,K)=ARRAY (1, J) 
ARRAY (L, J)=-SAVE 

DO 70 I=1,NORD 

IF (I-K)63,70,63 

ARRAY (I,K) =-ARRAY (1, K) /AMAX 
CONTINUE 

DO 80 I=1,NORD 

DO 80 J=1,NORD 

IF (I-K)74,80,74 

IF (J-K)75,80,75 | 
ARRAY (I,J) =ARRAY (1, J)+ARRAY (I,K) *ARRAY (K, J) 
CONTINUE 

DO 90 J=1,NORD 

IF (J-K)83,90,83 

ARRAY (K, J) =ARRAY (K, J) /AMAX 
CONTINUE 

ARRAY (K, K)=1. /AMAX 
DET=DET*AMAX 

DO 130 L=1,NORD 

K=NORD-L+1 

= ERC) 

TR G=K) 1d Ll 1.05 

DO 110 I=1,NORD 

SAVE=ARRAY (I,K) 

ARRAY (L,K)=-ARRAY (I,J) 
ARRAY (1,J)=SAVE 
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I=JK(K) 

IF (I-K) 130,130,113 

DO 120 J=1,NORD 
SAVE=ARRAY (K, J) 

ARRAY (K, J) =-ARRAY (I, J) 
ARRAY (I,J) =SAVE 
CONTINUE 

RETURN 

END 


QUADRATIC POLYNOMIAL FUNCTION 


FUNCTION QPOL(X,1,A) 
DIMENSION X(25),A(3) 

XI=X(1) 

QPOL=A (1)+A (2) *XI+A (3) *X1I**2 
RETURN 

END 


DERIVATIVE OF QUADRATIC POLYNOMIAL SUBROUTINE 


SUBROUTINE FDRIV(X,1,A,DLTA,NTRMS,DRIV) 
DIMENSION X(25) ,A(3),DLTA(3) ,DRIV(3) 


XI=X (I) 
DRIV(1)=1. 
DRIV(2)=XI 
DRIV(3)=XI**2 
RETURN 

END 


TOTAL FLUORESCENCE FUNCTION 


FUNCTION YTT1(A,T) 
DIMENSION A(3) 
F=A(1) 

F1=A(2) 

F2=A (3) 
YTTI=F+F1*T+E2*T**2 
RETURN 

END 
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BACKGROUND FLUORESCENCE FUNCTION 


FUNCTION YET1(A,DK,T) 
DIMENSION A(3) 

F=A(1) 

F1=A(2) 

F2=A (3) 
F3=(F1/DK-2.*F2/DK**2) *EXP (-DK&T) 
FG=F+F1*T+F2*T*%*2 
F5=-(F1+2.*F2*T) /DK 
F6=2.*F2/DK**2 
YET1=F3+F4+F5+F 6 
RETURN 

END 


DIFFUSION COEFFICIENT SUBROUTINE 


SUBROUTINE CLCD1(A,F,T,D) 
DIMENSION A(3) 
D=.1 
V=YET1(A,D,T) 
EF(E-V)1,4,3 
D=D-.01 

GO TO 9 
D=D+.0001 
V=YET1(A,D,T) 

IF (F-V)4,43 
WRITE(1,10)D,F,V 


FORMAT Cl) 0) (D=! E132 5a eho boa eo 


RETURN 
END 


v="",F8. 3) 
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